SACLAY-T09/158 
SISSA-64/2009/EP 



b 



Flavour violation in supersymmetric SO(10) 
unification with a type II seesaw mechanism 



Lorenzo Calibbi a , Michele Frigerio 6 , Stephane Lavignac 6 
and Andrea Romanino a 

a SIS S A/IS AS and INFN, 1-34013 Trieste, Italy 
Institut de Physique Theorique, CEA-Saclay, F-91191 Gif-sur-Yvette Cedex, France 



Abstract 

We study flavour violation in a supersymmetric SO(10) implementation of the type 
II seesaw mechanism, which provides a predictive realization of triplet leptogenesis. 
The experimental upper bounds on lepton flavour violating processes have a sig- 
nificant impact on the leptogenesis dynamics, in particular they exclude the strong 
washout regime. Requiring successful leptogenesis then constrains the otherwise 
largely unknown overall size of flavour-violating observables, thus yielding testable 
predictions. In particular, the branching ratio for fj, — ► e7 lies within the reach of 
the MEG experiment if the superpartner spectrum is accessible at the LHC, and 
the supersymmetric contribution to ek can account for a significant part of the 
experimental value. We show that this scenario can be realized in a consistent 
SO(10) model achieving gauge symmetry breaking and doublet-triplet splitting in 
agreement with the proton decay bounds, improving on the MSSM prediction for 
oi%{mz), and reproducing the measured quark and lepton masses. 
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1 Introduction 



Neutrino masses presumably arise at a scale much larger than the electroweak scale. If this is 
the case, a model-independent effective description of neutrino masses is possible at lower scales 
in terms of the dimension 5 operator hljh u h u / A [I], where l{ is the i-th. family lepton doublet, 
h u is the Y = +1/2 Higgs doublet and A is the scale at which the operator is generated, which 
can be as large as 10 15 GeV. While this effective description is essentially unique, the high- 
energy mechanism leading to the above dimension 5 operator is not. Generally, it is assumed 
to arise from the tree- level exchange of SM singlet fermions (type I seesaw mechanism [2]). In 
this case, the low-energy information from lepton masses and mixing only determines 9 out 
of the 18 high-energy seesaw parameters. Due to this arbitrariness, it is not possible to make 
definite predictions for the lepton asymmetry generated in the decays of the heavy singlet 
neutrinos [3], nor, in supersymmetric models, for the lepton flavour violating (LFV) effects 
induced by their Yukawa interactions [3]. 

On the other hand, the exchange of singlet fermions is not the only possible origin of 
the liljh u h u /A operator. In this paper, we consider the exchange of an SU(2)l triplet scalar 
with Y = ±1 (type II or triplet seesaw mechanism [5]). More precisely, we consider the 
SO(10) [6] implementation of the triplet seesaw mechanism proposed in Ref. [7J, in which all 
flavour parameters contributing to low-energy observables are determined in terms of the SM 
fermion masses and mixing^]. This allows to make testable predictions for these observables, 
up to a few unknown flavour-blind parameters. Particularly interesting is the possibility [7] of 
accounting for the baryon asymmetry of the universe via triplet leptogenesis. In usual type 
II models, it is necessary to introduce additional triplets or singlets in order to provide a rich 
enough flavour structure to induce a non- vanishing CP asymmetry in triplet decays [8]. This 
brings back into the game a number of flavour unknowns. On the contrary, in the scenario that 
we consider, the additional states are heavy quarks and leptons whose masses and couplings 
are determined in terms of low-energy parameters through SO(10) relations. The generated 
baryon asymmetry then directly depends on the light neutrino parameters. 

In this paper, we pursue the exploration of this scenario by studying the flavour- and 
CP-violating effects induced by the couplings of the heavy states to the MSSM squarks and 
sleptons. Assuming flavour-universal soft supersymmetry breaking terms at the GUT scale, 
flavour-violating observables are predicted up to a few unknown scale parameters and to a mild 
model-dependent uncertainty, thus allowing to test the scenario. Besides the contributions of 
the type II seesaw triplet and of its SU(5) partners already studied in Ref. [9], the presence of 
heavy quarks and leptons gives rise to additional contributions to the slepton and squark soft 
terms. In particular, they induce flavour and CP violation in the slepton singlet and squark 
doublet sectors, which were absent in the SU(5) type II seesaw model. 

The paper is organized as follows. In Section [21 we present the main features of the SO(10) 
scenario we consider. Section [3] contains a qualitative discussion of its predictions for flavour 
and CP violation. Model-building aspects including gauge coupling unification, doublet-triplet 
splitting and proton decay are briefly discussed in Section [H and addressed in greater detail in 
Appendices El E] and EH where the main ingredients of a realistic model are given. Section [5] 
presents our numerical results. Finally, we give our conclusions in Section[6l The superpotential 

1 Actually, the higher-dimensional operators needed to account for the measured masses of down quarks and 
charged leptons may affect this relation between high-energy and low-energy flavour parameters. We argue in 
Section [3]3] that the impact of this on physical observables is generally small, and neglect it in the following. 
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of the model, the boundary conditions for Yukawa couplings and soft terms and a subset of 
the renormalization group equations are displayed in Appendix iDl 

2 SO (10) unification with type II seesaw mechanism 

We consider a super symmetric SO(10) scenario with matter fields in 16 and 10 representations 
and type II realization of the seesaw mechanism, in which soft supersymmetry breaking terms 
arise at the GUT scale or higher. Since we are interested in the flavour- violating effects 
induced by the physics responsible for neutrino masses, we assume that the soft terms are 
flavour universal at the GUT scale. 

Let us first describe the field content. In terms of SU(5) multiplets, the three families of SM 
fermions are described by 5? © 10| , i = 1,2,3. In conventional SO(10) unification [6], 
5^ u<5) and 10^ u(5> are unified in a single 16^ O(10) together with a singlet (right-handed) neutrino 
participating in the generation of light neutrino masses through the type I seesaw mechanism. 
In this paper, we consider the alternative possibility that 10^ u(5> is embedded in a 16^ o(10) , 
while 5^ u(5) belongs to a 10^ O(10) (for similar or related approaches, see Ref. [H])- Having done 
this choice, we drop from now on the "SO(10)" superscript on SO(10) representations, and 
indicate the embedding of SU(5) representations into SO(10) representations by a superscript: 
for instance, 5*° is contained in 10j, while 5* 6 belongs to 16j. Besides the three 16j and 10j 
matter multiplets, the model also involves a 10, a 16 and a 54 Higgs multiplets. These are 
needed in order to generate the quark and charged lepton masses, as well as the neutrino masses 
through the type II seesaw mechanism. The corresponding SO(10) superpotential reads: 

W = Wvukawa + W seesaw + W G UT + W n onren. , (1) 
^Yukawa = ^ yy 1^16^10 + ^16^16 , (2) 

^seesaw = \ / ij 10 i 10 J 54 + 1 <t10 10 54 + i M 54 54 2 , (3) 

where WguTj whose explicit form is given in Appendices O and [Bj contains the terms respon- 
sible for SO(10) symmetry breaking and for the doublet-triplet splitting, and W nonren . includes 
the non-renormalizable operators needed to account for the measured ratios of down quark and 
charged lepton masses. The role of the 10, 16 and 54 representations is the following. The 
10 contains the Y = +1/2 MSSM Higgs doublet h u responsible for up quark masses. The 16 
plays a double role: it shares the Y = — 1/2 MSSM Higgs doublet with the 10 and (together 
with its 16 companion) it reduces the rank of the unified gauge group down to 4 through the 
vev of its SU(5)-singlet component; moreover, this vev pairs up the spare 5j 6 and 5j° and 
gives them a large mass, leaving only the 5*° massless. The 54 contains the SU(2)l triplet 
mediating the type II seesaw mechanism. The three singlet neutrinos contained in the 16j, on 
the other hand, do not couple to the light neutrinos at the renormalizable level. Assuming for 
definiteness that they acquire GUT-scale masses (e.g. from couplings 1^16^16 to three SO(10) 
singlets lj), we can completely neglect their contributions to the light neutrino masses, if any, 
and we are left with a pure type II seesaw mechanism. 

In writing Eqs. ([2]) and ©, we assumed that the superpotential is invariant under a matter 
parity, with the matter fields being the 16 j and the 10j. We also required the absence of mass 
terms of the form lOjlOj and of additional interactions that would induce a vev for the 54, 
in order to prevent a mixing between the 5*° and the 5* 6 . The absence of such a mixing 
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ensures that the light neutrino masses arise from a pure type II seesaw mechanism, and is a 
crucial feature of the scenario. The scale of light neutrino masses requires M54 to be smaller 
than the GUT scale; it is therefore natural to assume that this mass term arises at the non- 
renormalizable level, while it is forbidden or small at the renormalizable level. Furthermore, 
the leptogenesis scenario proposed in Ref. [7] requires that the 24 54 component of the 54 be 
heavier than the (15 © 15) 54 one. We refer the reader to Appendix [C] for the discussion of the 
splitting of the 54. 

The degrees of freedom surviving below the GUT scale are the MSSM states, which are 
massless before electroweak symmetry breaking, three heavy pairs of vector-like matter Ej 6 © 
5*°, and the components of the 54, which are also heavy. The light MSSM fields and the 
heavy ones are embedded in the SO(10) representations of the model as follows: 



(5) 



where dots stands for heavy fields which have been integrated out at the GUT scale (this 
includes in particular the coloured Higgs triplets mediating D = 5 proton decay). We have 
introduced an angle 6h parametrizing the mixing among Y = —1/2 Higgs doublets (0 < 9h < 
ir/2 with no loss of generality) and assumed that the Y = +1/2 MSSM Higgs doublet entirely 
resides in the 10, so that all components of the 16 are heavy (see Appendix [B] for an explicit 
realization of this). The SU(3)c x SU(2)l x U(l)y quantum numbers of the 54 components 
are A = (l,3)+i, £ = (6,l)_ 2/3 , Z = (3,2) +1/6 , S = (1,1) , T = (1,3) , O = (8, 1) and 
V ~ (3, 2)_ 5 / 6 , while the heavy L{ and Df obviously carry the same quantum numbers as the 
light k and d\. After breaking of the SO(10) gauge symmetry, the heavy quark and lepton 
fields acquire Dirac masses (Mi)c)ijDfDj + (ML)ijLiLj, where, at the renormalizable level and 
neglecting the renormalization group running below the GUT scale: 

(M2jc) - = (M L )ij = hijVi , (7) 

in which V\ is the vev of the SU(5)-singlet component of the 16. As for the SM fermion masses, 
they are given by (again neglecting corrections from non-renormalizable operators and the RG 
running of the Yukawa couplings): 

(M u )ij = yijV u , 

(M e )ij = sm6 H hijV d , (8) 
(M d )ij = smenhfjVd , 

and the neutrino masses are given by the type II seesaw formula: 

(M^ = -^/„. (9) 

Note that the down quark and charged lepton masses, which satisfy the SU(5) relation M d = 
Mj at the renormalizable level, are proportional to siylOh, the Y = —1/2 Higgs mixing 
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parameteiQ Due to the way the SM fermions are embedded into SO(10) representations, the 
up quark mass matrix is not correlated with the down quark and charged lepton mass matrices 
as in conventional SO(10) models, which accounts in a natural way for the stronger mass 
hierarchy observed in the up quark sector. 

The role of the SO(10) symmetry is to relate the masses and couplings of the heavy matter 
fields (Li, Li) and (Df,Df) to the ones of the light fermions, thus allowing to predict their 
contributions to observables such as the baryon asymmetry of the universe and the rates of 
flavour-violating processes. Indeed, Eqs. ([2])-([3]) and ©-([9]) show that all flavour param- 
eters involved in the superpotential, as well as the masses of the heavy matter fields, are 
determined by the SM fermion masses and mixings, up to flavour-blind factors and to the 
non-renormalizable contributions needed to make 7^ Mj. On the contrary, in conventional 
SO(10) models where neutrino masses arise from the type I seesaw mechanism, the Lagrangian 
below the GUT scale also depends on the flavour parameters encoded in the so-called R ma- 
trix [TT]. The SO(10) scenario studied in this paper therefore has a higher predictive power. As 
for the unknowns associated with the non-renormalizable operators in Wnonren.; we Wi h argue 
in Section 13.31 that they are unlikely to affect our results in a sizable way. 

3 Flavour and CP violation 

3.1 General structure of radiative corrections to the MSSM soft terms 

The /jjlOjlOj'54 interactions introduce a new flavour structure, directly related to the light 
neutrino masses, on top of the MSSM one. As noted in Ref. [12j . this can be considered as 
a truly minimal extension to the lepton sector of the minimal flavour violation hypothesis in 
the quark sector, which does not rely on flavour basis dependent assumptions. In any case, 
the /jjl0jl0j54 interactions give rise to new flavour- and CP-violating effects at low energy 
because of their well-known impact on the MSSM soft terms through radiative corrections [9]. 
In addition, the 16; 16,, 10 interactions induce new flavour and CP violation in the slepton 
singlet sector due to the presence of heavy lepton doublets in the 16j, as we are going to see. 

In order to compute these effects, we must write the renormalisation group equations 
(RGEs) for the MSSM parameters and integrate them. Before doing so, let us note that the Zj 
and Li fields can mix, since they have the same quantum numbers after breaking of the SO(10) 
symmetry (and similarly for the d\ and Df fields). Hence, the superpotential mass term for 
the lepton doublets reads, in full generality, Li(M£)ijLj + li(Mn)ijLj. The second term was 
omitted previously since Mil vanishes at the tree level, while Ml is given by Eq. ([7]). At one 
loop, however, wave function renormalization induces a small U-Lj mixing: 

(Afc)„ . (ML y^{^fi) . (10) 

This in turn necessitates a redefinition of the heavy and light lepton doublets, which affects 
all couplings involving lepton doublets. However, the size of the effect is small (especially for 
small or moderate tan/3), and we shall neglect it in the following. 

Before solving numerically the full 1-loop RGEs, let us illustrate the main features of the 
results by using the leading-log approximation for the exchange of the heavy degrees of freedom. 

2 This offers the possibility of explaining the smallness of the bottom/top mass hierarchy for moderate values 
of tan/3 in terms of a small Higgs mixing angle 9h ■ For reasons explained in Section |TT] however, we shall not 
use this possibility here. 
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In order to be able to identify each contribution, we assume that all components of the 54, 
namely S, T, O and the vector-like pairs (V, V) (A, A), (E, S) and {Z, Z) have different masses 
(the mass Mz of the [Z, Z) vector-like pair should not be confused with the Z boson mass 
mz)- We also consider the possibility that Mpc ^ Ml, but ignore the effect of corrections 
to the mass relation M^ = Mj. Let us denote by m 2 6 and m 2 the universal soft terms for 



the three families of 16j and 10j matter fields at the GUT scale; by TO54, fnie H and mio H the 
soft terms for the 54, 16 and 10 SO(10) multiplets; and by ao the universal A-term (defined 
by A x = aox, where x is any superpotential trilinear coupling x). Using the RGEs given in 
Appendix [HI we obtain (in matrix form): 



mf 



m 2 d c 



(mf)MSSU ~ ( 2m io + m 54 + «o) f j 



1 M 2 1 M 2 
_ l n GUT , £ 1 1V1 GXJT 

2 M 2 : " " 



M 2 



A 2 

M 2 



_i_ " l n j/i GUT + A l n ^GUT + - In ^GUT 
' M 2 +MlM* L 20 M 2 +MlM* L 2 My + M]^ C M^ C 

1 r A/T 2 /if 2 

' ' 1 • 2 ln GUT + In GUT 



("40MSSM ~ 7^32 ( 2m 10 + m 54 + Oo) / 



M 2 



+ ln 



M GUT , j ln jguj 



1 



M 2 
GUT 



_ l^jr + J_ ln 

3 M2+M], C M DC 15 M 2 +M|, C M 



f , (Ha) 



(lib) 



In contrast to the well-known type I seesaw case, the flavour-violating corrections to mf are 
determined by the light neutrino mass matrix, with no ambiguity due to high-energy flavour 
parameters. Therefore, while the absolute rate of a given LFV process is not known, the 
correlations between different LFV channels are predicted with little uncertainty (at least if 
the contributions of the slepton singlet sector, to be discussed below, are subdominant). The 
first term in the squared brackets of Eq. (jllb ) is induced by the seesaw triplet interactions 
and is present in all type II seesaw models. The second term, as well as the first two terms 
in the squared brackets of Eq. (jllb ). corresponds to the contribution of the SU(5) partners of 
the triplet [9] . The additional terms in mf and m\ c are due to the presence of heavy leptons 
and quarks in the 16j, and are characteristic of the model studied in this paper. In the limit 
M L = Mflc , M A = M E = M z , M s = M T = M = M v , the correction to the MSSM evolution 
is the same for mf and (m 2 ic ) T , as dictated by the SU(5) invariance of the interactions from 
which they arise. 

Another difference from other type II seesaw models comes from the corrections to the 
MSSM running of m 2 c and m 2 : 



m e 



m n 



(mgc)MSSM 
( m g)MSSM - 



cos 2 Oh 
(4ir) 2 

cos 2 9h 
(4vr) 2 



(2m 2 6 + m\ d + ag) y 
(2m 2 6 + m\ d + ao) y ] 



M 2 
2 In GUT 



M L M L 



M 2 
1 1V1 GUT 



M DC M. 



(12a) 
(12b) 



cos 2 #ff mf 0w + sin 2 ^ m 2 6fl . These are controlled by the up quark Yukawa 



couplings (evolved at the high scale). While the corrections to m 2 do not represent a deviation 
from the minimal flavour violation structure of the MSSM radiative corrections, the corrections 
to m 2 c are similar to the ones induced above the GUT scale by the top quark Yukawa coupling 
in SU(5) models [13J, although their origin is different. 
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3.2 Analytic approximations for the mass insertions and CP violation 



In order to be able to quickly estimate the size of various flavour- and CP-violating observables, 
it is useful to provide analytic expressions for the mass insertions parameters (i ^ j) [14]: 

/re \ — ( m f)ij /re \ - ( m %)ij /re \ _ (A e )jjV d 

{d LL )ij = -=g- , [d RR ) i3 = — ^— , [b RL ) i3 = =r^r , ("J 

e L e R C L 

where mg £ ijne R ) is an average doublet (singlet) slepton mass, and analogous quantities are 
defined in the up and down squark sectors. These parameters can be straightforwardly derived 
from Eqs. (jlip and (|12p . but some care is needed regarding CP-violating phases. Even if we 
take real boundary conditions for the soft supersymmetry breaking parameters at MgijT) we 
will end up with complex soft terms at the weak scale because their RGEs involve complex 
couplings. These couplings contain, in addition to the CKM and PMNS phases, extra CP- 
violating phases inherited from the SO(10) structure. In order to identify the latter, it is useful 
to write the SO(10) Yukawa couplings in an appropriate basis for the SO(10) matter multiplets 
16 j and 10 j, namel}H: 

y = UjD y U q , h = D h , f = U*D f rf , (14) 

where 

U q = Diag( e i *i,e i *2,e l *3)yDiag(e i *Se i ^,e i *3) . (15) 

In Eqs. (fT4"|) and (fl"5j) . D y , D^, Dj are real diagonal matrices, U and V are the PMNS and 
CKM matrices in the standard PDG parametrization (extended to include two "Major ana" 
phases p and a in the PMNS case), and <J>", <&f (i = 1,2,3) are extra SO(10) phases, five of 
which are independent (one can impose e.g. $3 = 0). For simplicity, we neglect the effects 
of the non-renormalizable operators needed to correct the mass relation Mj = Mj, and we 
therefore also assume Mpe = Ml = D^Vi- Below the GUT scale, we are free to rephase 
independently the MSSM fields so that the Yukawa couplings \ u ,d,e and /a only contain CKM 
and PMNS phases; however, the extra SO (10) phases will reappear in other couplings, such 
as Xd and fz, which enter the RGEs for the MSSM soft terms (we refer to Appendix ID1 for 
the definition of the superpotential couplings below the GUT scale). The outcome of this is 
that the mass insertion parameters depend on the phase differences — and <f>f — <3?j in 
addition to the 4 low-energy phases 5qkm, &pmns-> P an d c. 

We are now ready to write the mass insertion parameters (^jv)ij [M, N = L,R) in the 
leading-log approximation, assuming for simplicity Ma = Mx = Mz = M15, M$ = Mt = 
Mo = My = M24, a common soft supersymmetry breaking mass mo for all SO(10) chiral 
multiplets and, as before, a common A-terms ao for all superpotential trilinear couplings. We 
obtain: 

l^\8l L ) l3 « _ 3r»g + |ao| 2 g / 



16*\5 RR ) 13 « _ 2(3mg + |a | 2 ) o< (+ g -+j ) „ 



{Ct-T , (17) 



l^\5 RL ) l3 « -=^r- m ei CU^-^(CtrTm ej , (18) 



2me L me R 



In the following discussion, we neglect the effect of radiative corrections in the Yukawa sector. This allows 
us to disentangle the CKM and PMNS phases from the extra SO(10) phases, which would otherwise be mixed 
by the renormalization group running. 
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+ e 
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(20) 
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3a 
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MSSMjiJ 
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(22) 

(23) 
(24) 



where the dependence on the low-energy flavour parameters and phases is encapsulated in the 

(C)VLbB:\l ,, , • ■ . ,, . (/ ~ ;y 



coefficients (C^ SM )y, (C M g gM )y, C/j, C^ u and -D^™. The first two correspond to the CKM- 



induced MSSM contributions and are given by: 



mssm)ij — A 2 ^ti^* 



« ln (^Sf ) > (^sm), = £ fc A 2 fc y, fc ^ In , (25) 



where terms suppressed by \ c j\ have been neglected in (Cmssm)^- The other three coefficients 
correspond to the contributions of the heavy states present below the GUT scale and also 
contain a dependence on their masses. The one that involves the seesaw couplings is given by: 



rJ 



M GUT 



+ 5 faif< 



1 / Mq UT 



M GUT 



+ TT, k ,ifkflU ik U* l Y, a U ak U* al \n 



(26) 



where ft = (2M A /<rt> 2 ) and the Ms a 



1,2,3) are the common masses of the heavy 



quarks and leptons, and the coefficients that involve the up quark Yukawa couplings read: 



* C os 2 ^ A 2 ^ ln(^§F 



D 



i.i 



A£2 

M'i 



(27) 



where terms suppressed by X c /Xt have been neglected in C^ u . Due to this (very good) approx- 



imation, only the CKM and PMNS phases appear in Eqs. (|25p to (1271) . 

Let us have a closer look at the mass insertion parameters (|18|) to (|24[) , In the up squark 
sector, the new contributions are subdominant with respect to the MSSM corrections, which 
are enhanced by tan 2 /? and by a large logarithm. The situation is more interesting in the down 
squark and charged slepton sectors. In the former, the new corrections controlled by the top 
quark Yukawa coupling are again subdominant with respect to the MSSM corrections, while 
they lead to RR flavour violation in the latter. The new corrections induced by the seesaw 
couplings, on the other hand, are not suppressed by small CKM angles and can give rise to 
large flavour violations in the RR down squark sector and in the LL slepton sector [9]. One 
observes the following correlations: 
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d R 
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(4r) 

(Shi) 



'J 



rnlML) 



* i(<5> d -<5> d ) 
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2Rm\{8 d LL ) 



,■ e 1 J 



(28) 
(29) 
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5.8 x 10~ 5 
4.3 x 10" 4 


2.1 x 10~ e 
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1.4 x 10" 4 
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4.6 x 10~ 7 
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Table 1: Numerical values of the leptonic mass iusertiou parameters , ^b^r for two different 
values of sin 2 #13 and tan/3 and for the choice of parameters mentioned in the text. 



where R = cos 2 H ln(M^ UT /Mf 3 )/ [ ln(M^ UT /m|) + cos 2 ^ ln(M^, UT /M 5 2 3 ) ] . The first one 
is characteristic of the SU(5) extension of the type II seesaw mechanism, while the second one 
is analogous (although from a different origin) to the one arising from the running between 
Mp and Mqtjt in the minimal SU(5) model. In the case that the seesaw-induced corrections 
dominate, one further has: 

(Sunk « {Sll)h , ^ll)h « « R h . (30) 

m d R m d L ( 6 RL)ij = "Jg^meiO&k)*^** - ^ ■ ( 31 ) 

Eqs. (|28p . (|29p and (|3ip correlate the size of flavour-changing neutral currents (FCNCs) in 
the lepton and B/K sectors, as well as the charged lepton and quark electric dipole moments 
(EDMs), as we are going to see. 

For illustration, we give in Table [1] numerical values for the leptonic 6's. These values were 
obtained by numerically solving the RGEs for the following choice of high-energy parameters, 
motivated by the leptogenesis scenario of Ref. 0: Ma = Ms = Mz = M15 = 10 12 GeV, 
M s = M T = M = My = M 24 = 10 13 GeV, a = 3.4 x 10~ 2 , Vi = M GU t = 2 x 10 16 GeV and 
tan#fl- = 1. For the heavy (5j,5j) masses, we used M^ = m ei Vi/ sinOjiVd, neglecting SU(5)- 
breaking corrections. The supersymmetric parameters were chosen to be mo = M|/ 2 = 1.5 TeV, 
ao = 0, and tan/3 = 10 (50), yielding the following superpartner masses (for tan/3 = 10): 
Mw) = 112 GeV, M-+ = 218 GeV, M- Q = 723 GeV, sleptons around (1.5-1.6) TeV and first two 
generation squarks around 1.8 TeV (as well as fj, = 562 GeV). Finally, we took the best fit values 
of Ref. [15] for the measured neutrino oscillation parameters, together with mi = 0.005 eV and 
sin 2 #13 = (0.05) for the yet unknown parameters. All phases were set to zero. The resulting 
values of the leptonic mass insertions shown in Table [T] can be compared with the bounds 
coming from the non-observation of LFV decays of charged leptons (see e.g. Refs. |16l I17|, [18] ) . 
Although the sleptons are heavy in this example, the relatively large values of the 5's lead to 
a branching ratio for fi — > e^f just below the experimental bound for tan /3 = 10, and above it 
for tan/3 = 50. Using Eqs. (f28j) . (j29|) and (f31"|) . one can also compare the figures in Table [Q 
with the constraints on hadronic 5's coming from B and K physics (see e.g. Refs. |19} I18|). 

From Table [1] we can see that the LL mass insertions are the dominant source of lepton 
flavour violation in our numerical example, i.e. we are in the situation described by Eq. (|30p . 
This is due to the fact that the choice made for the model parameters leads to rather large 
/ij's (namely ^33 ~ 0.1): a smaller value of the ratio M^/a would yield smaller fij couplings, 
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hence smaller LL leptonic mass insertions, without affecting the RR mass insertions. We 
remind the reader that, while the overall size of the fy's is controlled by M/^/a, their flavour 
structure, up to small RGE effects, solely depends on the light neutrino parameters. Hence, 
even though the predictions for LFV processes span many orders of magnitude due to their 
strong dependence on the unkown scale parameter M^/a, the ratios of rates for different 
flavour channels are predicted with a mild dependence on the high-energy parameters, as long 
as the corrections controlled by the seesaw couplings dominate. This is precisely the case in 
our numerical example, in which the dominance of the LL mass insertions allows to estimate 
the ratios of LFV decays, as in usual type II seesaw models [§J: 



(32) 
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(33) 



0.17. 

The present upper bound on BR(n — > ery) thus excludes the possibility of observing r — > e^y 
in foreseeable experiments, while r — > ^7 may be accessible at a super-B factory [20] if ^13 is 
small. These ratios may change if the SO(10) parameters controlling the masses of the heavy 
states are varied but, as long as Eq. (|30|) is satisfied, the dependence is only logarithmic, while 
the branching ratios themselves are roughly proportional to (Ma/<j) 4 . 

Let us now turn our attention to CP violation. It is well known that, even if the "flavour- 
independent" phases carried by the diagonal A-terms and by the mu term (in the convention 
that the universal gaugino mass parameter is real) are absent, significant contributions to 
the leptonic and hadronic EDMs may arise from phases carried by flavour non-diagonal soft 
terms \21 \ [T6 l [22] . In the mass insertion approach, this corresponds to multiple insertions of 
<5's. For down quark and charged lepton EDMs, the dominant "flavour-violating" contribution 
(coming respectively from gluino and bino diagrams) is proportional to: 



Im 



(*L)« (*L)33 (s f RR h 



-/U tan /? m/ 3 Im 



(f = d,e), (34) 



in which the suppression due to the double flavour change is compensated for by the tan /3 m f 3 
enhancement. As for the up quark EDM, since (S RR )ij ~ its flavour-violating contributions 
are suppressed by small Yukawa couplings or CKM angles, making it much smaller than the 
down quark EDM. Together with Eqs. (|28j) and (129j) . this implies a correlation between the 
electron and neutron EDMs, which depend on the same combination of phases and flavour 
couplings: 



Im 



An 



°13 U 13 



(35) 



Using formulae available in the literature (see e.g. Ref. [23]), one can estimate the size of the 
charged lepton and neutron EDMs for the choice of parameters made in Table HJ With the 
superpartner masses given above and an average slepton (down squark) mass mg = 1.56 TeV 
(m,- = 1.8 TeV), one obtains, for tan/3 = 10: 



d ei ~ (-1.3xl0- 24 ecm)lm[(<5£ L ) 43 (^ K )3,], 

d n ~ (1 ± 0.5) (-2.9 x 10" 22 e cm) Im[(^ L ) 13 (O31 



(36) 
(37) 
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For the neutron EDM, we used the formula d n = (1±0.5) [lA(d d -0.25d u )+l.le(d c d +0.5d c u )} [M], 
in which we only kept the dominant gluino contributions to down quark electric and chromo- 
electric dipole moments, while we neglected the up quark dipole moments. The values of 
the <5's given in TableU together with 1(^)13(^)31! = 2 - 2 x 10 ~ 7 ( for #13 = °) and 
3.5 x 1(T 6 (for sin 2 9 13 = 0.05), lead to d e < (7 X 10" 33 - 10~ 31 ) ecm, d M < 3 x 10" 30 ecm and 
d n < (1±0.5) (6 x 10~ 29 — 10~ 27 ) ecm, where the upper ranges correspond to varying #13 from 
to its experimental upper limit. These bounds are well below the present experimental upper 
bounds (respectively 1.6 x 10~ 27 ecm [25], 1.9 x 10~ 19 ecm [26] and 2.9 x 10 _26 ecm [27]), while 
some tension with the constraints from LFV processes is already present for the same values 
of the 5's. Therefore, taking into account the FCNC constraints, we do not expect significant 
flavour-violating contributions to EDMs in the scenario studied in this papei0. 

A more promising CP-violating observable is the indirect CP violation parameter in the 
kaon sector, ek- A recent evaluation of the Standard Model prediction for ek suggests that 
its measured value still allows for a significant supersymmetric contribution [28J. In the mass 
insertion language, the relevant contributions are Im.[(5f IN )i2(SpQ)i2], with (M,N),(P,Q) E 
{(L, L), (R, R)} or {(L, R), (R, L)}. Since the LR mass insertions are suppressed by the down 
and strange quark masses and | {Sf L )i2\ is suppressed by | V^V^I, the leading contributions to ek 
(for / 33 > 0.01 - 0.03) are proportional to Im [{8 d LL )i 2 {5 d RR ) 12 } oc R' x Im [e***?-**) (C( 2 )* C&] 
and to lm[(5 RR ) 2 12 ] oc Im [e 2 ^*? - **) (C/2) 2 ], which depend both on CKM/PMNS phases and 
on the unknown extra SO(10) phases. As discussed at the end of Section[5J this might account 
for a significant part of the measured value of ek- 



3.3 Model dependence related to higher-dimensional operators 

Throughout the above discussion, we assumed that Eqs. ([7]) and (jHJ) hold at the GUT scale, 
yielding the following relations between the masses of the heavy (L,,Zj) and (Df,Df) fields 
and of the SM quarks and leptons: 

M L = . ^ M e , M DC = . ^ Mj , (38) 
sm 6 H v d sin 6 H v d 

as well as M d = Mj. The last relation is known to be in gross contradiction with the measured 
masses of the first two generation fermions, but since the corresponding Yukawa couplings are 
small, this can easily be cured by higher-dimensional operators involving fields with SU(5)- 
breaking vevs. These operators also affect, in a model-dependent way, the relations (|38D 
between the masses of the heavy and light matter fields. Moreover, they generally introduce a 
mismatch between the SO (10) bases in which M d and M e are diagonal, as well as between the 
SO(10) bases in which M d and Mqc (M e and M£) are diagonal. All these effects have an impact 
on the running of the MSSM soft terms and on the resulting flavour violation. However we 
argue below that, under reasonable assumptions, their impact on physical observables should 
remain small. 

The first effect, i.e. the fact that the masses of the (Lj, Lj) and (Df, Df) fields are no longer 
determined by the light fermion masses, introduces only a mild model dependence. Indeed, the 
heavy masses and Md? enter Eqs. (fTT|) and ([T2]) [or equivalently Eqs. ([26]) and ([27]) ] only 
logarithmically. On top of that, in the case that M24 > Ml 2: dc, mf and m dc only depend on 



4 Of course, "flavour-independent" phases carried by diagonal A-terms or by the mu term may still induce 
large EDMs. 
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Ml 3 and Mpc , which are presumably little affected by the higher-dimensional operators, while 
in the case that M24 > Ml\,d5 some model dependence arises from the contributions of these 
operators to M^ 2 and Mpc. As for m%, it is dominated by M^-dependent contributions. 

The second effect, i.e. the mismatch between the mass eigenstate basis of the heavy and 
light quark (lepton) fields, introduces new flavour parameters in the model in the form of 
unitary matrices connecting a given SO(10) basis to another. As a result, the couplings of 
the components of the 10j (namely the Zj, d% Li and Df fields) to the components of the 
54 are no longer equal at the GUT scale, when expressed in terms of mass eigenstates. For 
instance, /a and fz (see Appendix [D] for the definition of these couplings) now differ by a 
(model-dependent) unitary matrix. It follows that the radiative corrections to the MSSM 
soft terms do not only depend on low-energy flavour parameters, but are also sensitive to 
these high-energy flavour parameters. However, since the unitary matrices originate from the 
higher-dimensional operators correcting the masses of the first two generation down quark and 
charged lepton, it is natural to expect that they are characterized by small mixing angles, 
thus only mildly affecting the flavour structure of the couplings, which on the contrary is 
characterized by large 612 and #23 angles. We shall therefore neglect this model dependence in 
the following, and simply input the SU (5) relations (|38p in our numerical studjd. 

For completeness, we list below the D=5 operators which can correct the mass relation 
M d = Mj, assuming the field content and pattern of vevs of Appendix lAl 

(i,2) (16il6)io(10j45i,2)io ^(i,2) (16 i 45i, 2 )i6(10jl6)ifl (16a6)io(10 J 54 / ) 10 

a ij a ' *i A ' Cij A ' [ ' 



where the subscripts specify the contraction of SO(10) indices and A is the cutoff. In Eq. 
the non-vanishing vevs of 45i and 45 2 are along the T^r and B — L directions, respectively, 



l(2) 

J ij ! 

the following corrections to the mass matrices of the heavy and light matter fields: 



and 54' has a vev in the Pati-Salam singlet direction. The operator b\- , for instance, yields 



(M^j = sin 9 H v d (h^ -3cb^y (M L ) i:j = V 1 (h ij +3cbW 

(A#)ij = sm6 H v d (h^ - cbf)) , (M D c) tj = V l (h^ + cbf) 



(40) 



where c = 2^2/3 Vb-l/A. Even in the simple case where a single operator (M 2 ) in this 

example) is present, the experimental values of the down quark and charged lepton masses 

(2) 

are not sufficient to fully determine the couplings hij and b- . Hence, the heavy field masses 
and the mixing patterns in the four mass matrices are model dependent. As argued above, 

however, this is unlikely to affect the radiative corrections to the MSSM soft terms in a sizable 

(2) 

way if the couplings hij and b\j have a hierarchical flavour structure. 
4 Model-building aspects 

Before presenting our numerical results for flavour violation in the next section, let us briefly 
present the main ingredients necessary to promote the SO(10) scenario studied in this paper to 
a realistic model. This includes the dynamics of gauge symmetry breaking, a viable doublet- 
triplet splitting mechanism and the generation of the intermediate scales associated with the 



5 We note in passing that there is enough freedom in the higher-dimensional couplings (see the example 
below) to account for the measured down quark and charged lepton masses without introducing new flavour 
parameters in the model, i.e. with Md, M e , Md<= and Ml all diagonal in the same SO(10) basis. 
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components of the 54 multiplet consistently with gauge coupling unification. More details can 
be found in Appendices EJ [B] and [Cj 

4.1 Gauge symmetry breaking and doublet-triplet splitting 

In order to break the SO(10) gauge group down to SU(3)c x SU(2)l x U(1)y, we introduce 
two adjoint Higgs representations 45 1 and 452 with non- vanishing vevs in the T^r and B — L 
directions, respectively, and a 54' which breaks SO(10) down to its Pati-Salam subgroup (we 
also introduce an SO(10) singlet). The rank of the gauge group is broken by the (16, 16) pair. 
A superpotential that enforces this pattern of vevs is given in Appendix(Aj In this example, all 
vevs are proportional to a single mass parameter, Miq. Hence, if all superpotential couplings 
are of order one, SO(10) is broken down to SU(3)c x SU(2)l x U(1)y in one step. Note that 
the 54' cannot be identified with the 54 multiplet involved in the seesaw mechanism, since 
a vev of the latter would induce an unwanted 5f°/5} 6 mixing (see Section [2]), and also for 
reasons related to the doublet-triplet splitting mechanism. 

In order to avoid proton decay at an unacceptable rate, we must split the electroweak 
doublets from the colour triplets in the 10 and 16 Higgs multiplets. This can be achieved by 
a generalization of the missing vev mechanism [29] involving an additional 10' Higgs multiplet 
and an adjoint Higgs representation with a vev aligned along the B — L direction, which is 
precisely the role of the 452 (see Appendix|B]for details). After SO(10) symmetry breaking, all 
colour triplets acquire GUT-scale masses, and a single pair of Higgs doublets (h u , h^) remains 
massless, the Y = +1/2 electroweak doublet in 10 and a combination of the Y = —1/2 
electroweak doublets in 10 and 16: 

h u = H u ° , h d = cos 9 H Hf + sinO H Hf , (41) 

where the Higgs mixing angle is given by tanflff = fjVi/Mie , with fj a superpotential coupling. 
As discussed in Section [21 sin Oh ^ is necessary to give masses to the down quarks and 
charged leptons, which live in the lOj's. Furthermore, the leptogenesis scenario of Ref. [7] 
assumes that contains a non- negligible component (cos 6 jj 7^ 0). 

Since the down quark and charged lepton masses are proportional to sin 9h at the renor- 
malizable level, the top/bottom mass hierarchy could in principle be explained by a small 
Higgs mixing, rather than by a hierarchy in the Yukawa couplings or a large value of tan/3. 
This, however, would imply rather heavy vector-like matter fields, since according to Eq. (138ft 
their masses are inversely proportional to sinOn- On the contrary, leptogenesis requires a rel- 
atively light (Li,Xi) pair in order to satisfy the condition for a non- vanishing CP asymmetry 
(Ma > 2Ml x ). We shall therefore stick to the case sinOn ~ 1. 

4.2 Proton decay 

Although all colour Higgs triplets acquire GUT-scale masses through the doublet-triplet split- 
ting mechanism, their contribution to the proton decay mode p — * K + V can still exceed the 
experimental limit. To suppress it further, one must impose additional constraints on the 
combination of doublet-triplet splitting parameters which enters the proton decay rate. More 
precisely, one must either arrange some cancellation among these parameters, or allow a pair of 
Higgs doublets to lie at an intermediate scale Mh < 10 14 GeV (see Appendix IB. 21 for details). 
In this paper, we choose the latter option, which in turn affects the running of the gauge cou- 
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plings and leads us to split the (15 © 15) 54 and 24 54 multiplets in order to restore successful 
gauge coupling unification, as we discuss below. 

4.3 Intermediate scales and gauge coupling unification 

The leptogenesis scenario of Ref. J7] requires the 15 © 15 component of the 54 to be lighter 
than its 24 component. Actually the full SU(5) multiplets are not needed for leptogenesis: 
only the (A, A) pair in 15 © 15, which is responsible for the type II seesaw mechanism, and 
the S and T fields in 24 (or just one of them) are really necessary. The main reason for using 
complete SU(5) representations in Ref. [7] was gauge coupling unification. However, we have 
seen that the experimental constraint on the p — ► K + v rate can be satisfied by allowing a pair 
of Higgs doublets to lie at an intermediate scale, which in turn spoils unification. As shown in 
Appendix[Cl this can be cured by splitting the masses of the components of (15 © 15) 54 and 
24 54 in an appropriate manner. This also has the advantage of maintaining perturbativity 
above the GUT scale, and of avoiding too strong flavour-violating effects in the region of 
the parameter space relevant for leptogenesis, where the couplings are large and induce 
important radiative corrections to the MSSM soft terms. 

Two simple possibilities emerge from the analysis of gauge coupling unification done in 
Appendix [C] (the ratio Mt/Ma = 10 is motivated by leptogenesis): 

• model (i): intermediate (A, A), (£, S) and T, with 

M A = M E , M T = 10M A , M H = 10 14 GeV , (42) 

• model (ii): intermediate (A, A), (£,£), S, T and O, with (S is actually irrelevant for 
gauge coupling unification) 

Ml/ 'Ml « 10 15 GeV , M s = M T = M Q = 10M A , M H = 10 14 GeV , (43) 



all other components of the 54 having GUT-scale masses. The desired splittings can be 
achieved by appropriate higher-dimensional operators (see Appendix [C]) . In model (i), the 
unification scale lies almost one order of magnitude lower than in the MSSM, which is at odds 
with the experimental limit on p — > 7r°e + . While this might be cured by 2-loop running and 
GUT threshold corrections, we prefer to adopt model (ii) for the numerical study of the next 
section. We note in passing that both models give similar (although quantitatively different) 
results for flavour violation. For illustration, let us consider model (ii) with M A = 10 12 GeV, 
M s = M s = M T = M = 10 13 GeV and M H = 10 14 GeV (the other parameters are chosen 
to be V\ = A/qut, Aff = |<r A (/u = Ma)| = 0.045, tan/3 = 10 and tan#fl- = 1). The spectrum 



of heavy states below the GUT scale is shown in Fig. 1(a), and the renormalization group 



running of the gauge couplings in Fig. 1(b) . At the 1-loop level, one finds: 



M GU T = 1-2 x 10 16 GeV , « G ut = 1/12.5 . (44) 

The prediction for a^(mz\ including super symmetric thresholds and the 2-loop MSSM run- 
ning, agrees with the measured value within la. Finally, the Landau pole lies one order of 
magnitude above Mqtjt- 
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Figure 1: (a) Spectrum of heavy states below the GUT scale in model (ii) for the following 
choice of parameters: M A = 10 12 GeV, M s = M s = M T = M = 10 13 GeV, M H = 10 14 GeV, 
V\ = MquTi = \ a A.(^ = = 0.045, tan/3 = 10 and tan Oh = 1. (b) Renormalization 

group running of aj x {ij) (i = 1, 2, 3) between p = 10 10 GeV and p = M G tjt- 



5 Numerical results 

In this section, we perform a numerical analysis of flavour violation in the SO(10) model defined 
above. Since the flavour-violating effects we want to study arise from radidative corrections to 
the sfermion soft terms, we numerically solve the full 1-loop RGEs from the unification scale 
Mqut down to the low-energy scale -Mewsb = ^ m t 1 m t 2 > taking into account the presence 
of additional states and interactions below the GUT scale, which modify the running of the 
MSSM soft terms as outlined in Section [3j A subset of these RGEs is given in Appendix |D| 
where the definitions of the various superpotential couplings and soft terms can also be found. 

In order to set the initial conditions at the GUT scale for the superpotential couplings and 
to compute the masses of the heavy (5j, 5j) pairs (which depend on the SO(10) couplings 
hij), we first evolve the Yukawa matrices and the low-energy neutrino mass operator from low 
energy up to the seesaw scale Ma, where the couplings are computed according to the 
seesaw formula (M u )ij = —XHfijVu/^M^, in which Xh = IcaCa 4 = -^a)|- Then we evolve the 
fij together with the other Yukawa couplings up to the GUT scale. This procedure is iterated 
until convergence of the heavy (5j, 5j) masses is reached. The other mass scales are taken 
as inputs and varied in order to study their impact on leptogenesis and on the flavour-violating 
observables listed in Table[2j Moreover, to satisfy the proton decay constraint consistently with 
gauge coupling unification, we split the components of the 54 as indicated in Eq. (l43j) . Thus, 
specifying the triplet mass Ma is enough to fix the masses of all 54 components. The value of 
Mh has a weaker impact on gauge coupling unification, and we keep it to be 10 14 GeV. The 
other inputs in the procedure, apart from the SM and supersymmetric parameters, are the 
seesaw coupling A#, V\ and tan0#, which are needed to compute the heavy (5j, 5j) masses. 

In order to single out the flavour- violating effects arising from radiative corrections, we 
assume universal boundary conditions for the soft supersymmetry breaking terms at the GUT 
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Observable 


Bound 


rtei. 


BR(fi — >■ cy) 


< 1.2 x 10~ n 
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BR(/i — > eee) 


< 1.0 x 1(T 12 
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CR(/x -► e in Ti) 


< 4.3 x 1(T 12 
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BR(r -> 


< 4.4 x 1(T 8 
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BR(r -» e7) 


< 3.3 x 10' 8 


33j 


£K 
AniK 
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Am Ss 
BR(6 -> 87) 
BR(£? S -> /x/x) 
Amg 
Am/) 


(2.229 ±0.012) x 10~ 3 

(3.483 ± 0.006) x 10~ 12 MeV 

(17.77 ± 0.12) ps 

2.77 x 10~ 4 - 4.33 x 10" 4 (3d) 

< 4.7 x 10~ 8 

(0.507 ± 0.005) ps" 1 
(0.0237^^6) ps -i 


[34j 

m 

10 A 1 

[34] 
[35] 
[36] 
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34j 



Table 2: Flavour-violating observables studied in the numerical analysis and their experimental 
values or upper bounds. 
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Figure 2: Feynman diagrams responsible (together with the CP-conjugated diagrams) for the 
CP asymmetry in the decays of the scalar triplet into Standard Model lepton doublets. 



scale, with common gaugino mass parameter Mj/2, scalar soft mass thq and ^4-term ao- For 
definiteness, we set ao = in the following. In a later stage, we shall also comment about the 
effect of relaxing the equality of the soft mass parameters for different SO(10) multiplets, still 
assuming flavour-blind supersymmetry breaking. After having performed the running from 
the GUT scale to low energy, we check that electroweak symmetry breaking does take place 
and that no tachyonic states are present in the spectrum. We then compute the masses of all 
Higgs bosons and superpartners and impose the mass limits coming from direct searches at 
LEP and at the Tevatron. For the rates of LFV processes, we use the expressions of Ref. [37]. 
The supersymmetric contributions to BR(B® S — ► /i + /x~) are estimated using the formulae of 
Ref. [38j, while BR(6 — ► sj) is computed with the help of the routine SusyBSG [39]. The 
supersymmetric contributions to the meson mass splittings Amx, Amp, Amg, Amg s and to 
the indirect CP violation parameter ek are computed in the mass insertion approximation, 
using the formulae of Ref. [30] ■ We recall in Table [2] the experimental values and upper limits 
for these observables. 

As observed in Section 13.21 and illustrated in Fig. [31 the predictions of the model for 
LFV processes span many orders of magnitude, due to their strong dependence on the ra- 
tio Ma/A//. In order to obtain definite predictions, we shall restrict the seesaw parame- 
ter space to the region favoured by the leptogenesis scenario of Ref. [7J, which is built in 
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1Q 11 1Q 12 1() 13 10 14 1() 15 

M A (GeV) 

Figure 3: Dependence of the CP asymmetry £a on Ma and Ah for tan f3 = 10, V\ = Mqut 
and maximal Higgs mixing (tan 9h = 1)- The unmeasured neutrino parameters are chosen to 
be mi = 0.005 eV, sin 2 #13 = 0.05, p = 7r/4, a = tt/2 and 5 = 0. The area below the white line 
does not satisfy the conditions for a large efficiency 

the SO(10) model considered here. The source of the lepton asymmetry is the CP asym- 
metry in triplet decays, which arises from the interference between the tree-level and one- 
loop diagrams shown in Fig. [2] and is given by (in the limit <C Ma < M^ 1 + M^ 2 ): 
e A ^ (0.1/107r)Im[/u(r//*)li] / [Tr(/D + |/n| 2 + I A^| 2 (l + cos 4 H )} , where M s = M T = 
IOMa has been used. The presence of heavy lepton fields with hierarchical masses is crucial 
for generating a non-vanishing CP asymmetry. In particular, the condition Ma > 2M^ must 
be satisfied for the loop integral to have an imaginary part. 

Let us identify the region of the (Ma, Ah) parameter space in which successful leptogenesis 
is possible. One can distinguish between two regimes [7J: the first one, characterized by a 
large value of the CP asymmetry in triplet decays and by a strong washout, requires order one 
values of the fij couplings. This however is in conflict with the experimental upper bounds 
on lepton flavour violation, unless the super symmetric spectrum is very heavy. In the second 
regime, the CP asymmetry is smaller, but the efficiency factor accounting for the dilution of the 
generated lepton asymmetry by washout processes can be of order one. Given that the observed 
baryon asymmetry is reproduced for rye a — 10~ 8 , where rj is the efficiency factor, successful 
leptogenesis is possible for £a > 2 x 10~ 8 . In Fig. El we plot £a as a function of the seesaw 
parameters Ma and Ah- The light neutrino parameters are chosen as in Ref. [7J, with the 
best fit values for the measured oscillation parameters taken from Ref. [15] and mi = 0.005 eV, 
sin 2 #13 = 0.05, p = vr/4, a = vr/2 and 5 = 0. In the region M A < 6 x 10 11 GeV, the 
condition for a non- vanishing CP asymmetry, Ma > 2M i ± , is not satisfied, while in the small 
Ah region some of the fij couplings become non-perturbative below the GUT scale, making 
the computation of £a non reliable. The white line separates the weak and strong washout 
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BR(n->e y) 




M A (GeV) 

Figure 4: Dependence of BR(/i — > cy) on Ma and Xh, for the following choice of mSUGRA 
parameters: ttiq = M\m = 700 GeV, ao = 0, tan (3 = 10 and /x > 0. The other parameters are 
chosen as in Fig. [3j The white dots indicate the area where both BR(/x — ► e'y) < 1.2 x 10~ n 
and e A >2x 10~ 8 . 



regimes: below this line, the conditions for a large efficiency, T(A — > L\L\) < H(M&) and 
T(A — ► H U H U ) > H(M/y), are not satisfied (see Ref. [7] for details). As shown by Fig. [3j £a 
can reach sizable values in the large efficiency region. 

Let us now see how the LFV decay [i — ► cy constrains the parameter space of Fig. [3j 
This is shown in Fig. U] for a point of the MSSM (mSUGRA) parameter space giving a rather 
light superpartner spectrum: ijiq = M\m = 700 GeV, ao = 0, tan/? = 10 and \x > 0. The 
brown (dark) area is excluded by the experimental upper limit BR(/j — > erf) < 1.2 x 10~ n . 
Comparing Fig. 2] with Fig. [3j one observes a tension between the requirement of successful 
leptogenesis and the experimental constraint on fx — > 67, due to the fact that both e& and 
BR(/i — > &y) grow with the ratio M^/Xjj (to which the couplings are proportional). Indeed, 
the region of the (Ma, Xh) parameter space consistent with both BR(/j — ► erf) < 1.2 x 10 -11 
and £a > 2 x 10~ 8 , marked with white dots in Fig. U is rather small for the chosen MSSM 
parameters^. An heavier superpartner spectrum and/or a smaller value of tan (3 would increase 
the size of this region by relaxing the /i — > e^f constraint, without affecting leptogenesis. A 
smaller value of $13 would also relax the \i — ► e7 constraint, but it would simultaneously reduce 
£A- Nevertheless, the requirement of successful leptogenesis (together with the non-observation 
of \i — ► erf) considerably restricts the seesaw parameter space, thus allowing us to make testable 
predictions for flavour-violating observables. 

In the following, we present our results for a set of seesaw parameters belonging to the region 

6 In passing, this shows that the supersymmetric version of the leptogenesis scenario proposed in Ref. [7] 
can be excluded on the basis of low-energy flavour physics measurements, if the supersymmetric partners are 
accessible at the LHC (and barring cancellation with other sources of flavour violation in the soft terms). 
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Figure 5: Predictions for BR(^ — ► &y) in the (mo, Mu^) plane, assuming ao = 0, tan/3 = 10 
and fx > 0. The parameters of the SO(10) model are: M A = 10 12 GeV, X H = 0.045, V\ = M GUT 
and tanOjj = 1. The neutrino parameters are chosen as in Fig. [3j 



where successful leptogenesis is possible, namely we take Ma = 10 12 GeV and Ah = 0.045, 
together with V\ = Mqut an d tan. Oh = 1 (which gives £a — 2.5 x 10~ 8 ). The corresponding 



spectrum of heavy states is shown in Fig. 1(a) As for the neutrino parameters, we choose 



them as in Fig. [3J We are now ready to present the predictions of the model for various 
flavour-violating observables as a function of the MSSM parameters. In Fig. [5l the contours 
of BR(/i — ► e*y) are plotted in the (mo, M]/ 2 ) plane for ao = 0, tan/3 = 10 and fj, > 
(the contours corresponding to different values of tan (5 can be easily deduced from Fig. [5] by 
noting that BR(/x — > cy) approximately scales as tan 2 /3). The parameter space is bounded 
by the LEP limits on superpartner and Higgs boson masses for low values of M1/2, and by 
the absence of radiative electroweak symmetry breaking for large values of mo and moderate 
Mi/2- BR(6 — > sj) gives a constraint similar to the Higgs mass bound, while the other hadronic 
observables of Table [2] do not significantly restrict the parameter space. This implies that the 
flavour physics signatures of the model are expected to show up in the lepton sector rather 
than in the hadronic sector (with the possible exception of ek discussed at the end of this 
section) . One can see from Fig. [5] that the on-going experiment MEG [41] , which aims at 
a sensitivity of 0(1O~ 13 ) on BR(/i — > e^y), will probe the model over a large portion of the 
MSSM parameter space. Interestingly, this region approximately corresponds to the one that 



is accessible at the LHC. This can be seen in Fig. 6(a) , where the contours of the gluino and of 



the lightest stop masses are plotted in the (mo,M 1 / 2 ) plane, and the area that will be probed 
by MEG is marked with dots. For completeness, we show in Fig. |6(b)] the full super symmetric 
spectrum corresponding to the mSUGRA point of Fig. [H i.e. mo = M 1 / 2 = 700 GeV, ao = 0, 
tan/3 = 10 and [i > 0. Note that this spectrum is characterized by rather heavy sfermions 
compared to neutralinos and charginos. This is due to the relatively large value of the unified 
gauge coupling, which enhances the gauge contributions to the running of sfermion masses at 
high energy. As a result, the lightest neutralino is found to be the LSP over the whole (mo, 
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Figure 6: (a) Contours of the gluino mass (blue dashed lines) and of the lightest stop mass (red 
solid lines) in the (mo, Mi/ 2) plane, for the same choice of parameters as in Fig. [5j The grey 
region is excluded by the constraints discussed in the text, including the experimental upper 
limit on BR(/i — ► ej). The dotted area will be probed by MEG. (b) Supersymmetric spectrum 
corresponding to the mSUGRA point uiq = Mi/ 2 = 700 GeV, ao = 0, tan/3 = 10 and fx > 0. 



M1/2) parameter space. For the spectrum showed in Fig. |6(b)[ the sleptons are heavier than all 
neutralinos and charginos and cannot be produced in cascade decays of squarks at the LHC. 
However, this possibility is recovered for Mi/ 2 3> tuq. 



Let us now present the predictions of the model for other LFV processes. In Fig. 7(a) 
BR(r — * n'y), BR(/j — * eee) and the fi — e conversion rate in the Titanium nucleus are plotted 
against BR(/i — ► cy) for the same parameter choice as in Fig.[5j BR(r — > ej) is not shown since, 
as discussed in Section [3T21 it is generally smaller than BR(r — ► ^7). Given the experimental 
upper limits shown in TableEJ BR(/i — ► ej) is at present the most constraining LFV observable. 
Since BR(r — ► /x7)/BR(/x —* erf) = 0(1), in agreement with Eq. (|32[) and with the correlation 
observed in type II seesaw models for large #13 [9], r — > /J7 is out of reach of super B factories, 
which are expected to achieve a sensivity of 10 -9 on its branching ratio |20j . fj, — e conversion 
looks more promising, given that proposed experiments at Fermilab [42] and at J-PARK |43j 
aim at respective sensitivities of 10~ 16 and 10~ 18 on CR(/i — > e in Ti). If approved, these 
experiments would test the model well beyond the MEG reach. 

In order to estimate the impact of non-universal (but flavour-blind) boundary conditions 
on the predictions for LFV observables, we performed a random scan of the soft scalar masses 
for different SO(10) multiplets between and 3TeV, assuming a fixed value M x / 2 = 700 GeV 
of the common gaugino mass parameter. The results for BR(/i — > ej) are shown in Fig. |7(b)| 
where the blue dotted lindll corresponds to the universal case mie = mio = rriiQ H = mio H = 
77154 = rriQ. Relaxing the universality of soft scalar masses can enhance or suppress BR(/i — * cy) 
by up to 2 orders of magnitude, but most of the points still remain within the reach of MEG 
(unless the lightest slepton is very heavy). 

As mentioned earlier, LFV observables provide much stronger constraints on the model 



7 In the universal case, radiative electroweak symmetry breaking does not take place for large mo values, 
which explains why the blue dotted line stops at mf 1 « 1750 GeV. 
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Figure 7: (a) BR(r — ► /ry), BR( / u — ► eee) and CR(/j, — ► e in Ti) versus BR(/j — ► &y) for the same 
choice of parameters as in Fig. [SJ scanning over m,Q and M 1 / 2 in the ranges < mo < 3 TeV and 
< M1/2 < ^ TeV. The grey area is excluded by the experimental upper limit on BR(/j — > e7). 
(b) BR(/j — ► e7) as a function of the lightest slepton mass for non-universal soft scalar masses 
Tn\e, mio, mie H , mio H and 77154 in the [0,3 TeV] range, and M 1 / 2 = 700 GeV. The blue dotted 
line corresponds to the universal case m\% = m\Q = m\§ H = mio H = 77754 = ijiq. 

studied in this paper than hadronic observables. A possible exception is represented by the 
indirect CP violation parameter in the kaon sector, ek- It is well known that ek is very sensitive 
to new sources of flavour and CP violation in the 1-2 down squark sector. In the present model, 
radiative corrections generate large off-diagonal entries both in the LL slepton mass matrix 
(which leads to large LFV rates) and in the RR down squark mass matrix. Moreover, order 
one phases in the PMNS matrix, hence in the couplings, are needed to account for the 
baryon asymmetry of the universe. Unfortunately, unknown extra SO(10) phases spoil the link 
between CP violation in the neutrino sector and the RG-induced CP violation in the sfermion 
mass matrices (see the discussion in Section f3. 2D . Nevertheless, barring cancellations between 
contributions carrying different phases, it is reasonable to expect a large imaginary part of 
(<5jyj)i2- In Fig. El the contours of the supersymmetric contribution to ek are plotted in the 
(mo, M1/2) plane for the same choice of parameters as in Fig. [5l assuming arg [(<5jyj)i2] = 0.5. 
While low values of M1/2 would give too large a contribution to a rather light superpartner 
spectrum can account for up to (10 — 20)% of its experimental value. According to Ref. [28J, 
this is precisely what is needed in order to reconcile the SM prediction for ek with experiment. 

6 Conclusions 

We have studied flavour violation in a supersymmetric SO(10) implementation of the type II 
seesaw mechanism, which provides a predictive realization of triplet leptogenesis. In this sce- 
nario, the high-energy flavour parameters involved in the computation of the bayon asymmetry 
of the universe and of flavour-violating observables are determined in terms of the Standard 
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Figure 8: Contours of |£i^ sy |/|£i^ P | in the (mo,Mi/2) plane (where e^ sy is the supersymmetric 
contribution to ek and e^ p its central experimental value), for the same choice of parameters 
as in Fig. [5] and assuming arg [(£^2)12] = 0.5. 

Model fermion masses and mixing, up to mild model-dependent uncertainties. The overall size 
of the FCNC effects is then controlled by a few unknown flavour-blind parameters, while the 
ratios of FCNC rates for different flavour channels, such as BR(t — > fj,^)/BR(fi — » ej), mainly 
depend on low-energy parameters. 

The features of flavour violation in the SO(10) scenario studied in this paper present some 
interesting differences with the SU(5) implementation of the type II seesaw mechanism [9], 
because of additional contributions coming from the heavy matter fields crucial for leptogenesis. 
These give rise, for example, to radiative corrections to the soft mass matrices m^ c and m 2 q 
controlled by the top quark Yukawa coupling. Moreover, the presence of a built-in leptogenesis 
mechanism provides a criterion for fixing the values of the unknown flavour-blind parameters, 
thus yielding testable predictions for the rates of flavour-violating processes. Imposing the 
conditions for successful leptogenesis together with the experimental constraints on FCNCs, 
we found that the predicted branching ratio for /i — > e7 lies within the sensitivity of the MEG 
experiment if the superpartner spectrum is accessible at the LHC, while r — > ^7 is out of reach 
of future super B factories, fj, — e conversion on Titanium is also a promising process, with 
a predicted rate within the reach of proposed experiments at Fermilab and J-PARK even in 
regions of the parameter space where BR(fi — ► e^y) lies below the MEG sensitivity. Hadronic 
observables only receive small contributions once the experimental bounds on LFV processes 
are imposed, with the possible exception of ek- 

We have also studied flavour- violating contributions to leptonic and hadronic EDMs, as well 
as to the ek parameter. The CKM and PMNS phases, together with new CP-violating phases 
associated with the SO (10) structure, can in principle induce sizable contributions to these 
observables, even if the soft terms are real at the high scale. The experimental bounds on LFV 
processes, however, prevent significant contributions to the EDMs, while a sizable contribution 
to ek (which might be necessary to account for its measured value) is still possible. 

The predictivity of the scenario studied in this paper relies on the SO(10) relations between 
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the flavour structures of the SM and heavy matter fields. These relations, however, only hold at 
the tree level and may be affected by model-dependent corrections from the non-renormalizable 
operators necessary to account for the measured quark and lepton masses. Nevertheless, the 
impact of these corrections on our results can be estimated to be mild. 

The presence of heavy states below the GUT scale, besides inducing low-energy flavour and 
CP violation through radiative corrections, also has an impact on the superpartner spectrum. 
Indeed, their contribution to the beta functions leads to a relatively large value of the unified 
gauge coupling, which enhances the gauge contribution to the running of sfermion masses. 
One of the consequences of this is that the lightest neutralino turns out to be the LSP in an 
unusually large portion of the parameter space. 

Finally, we have shown that the SO (10) implementation of the type II seesaw mechanism 
studied in this paper can be promoted to a consistent model including the dynamics of gauge 
symmetry breaking, doublet-triplet splitting consistent with the present experimental bounds 
on proton decay, and gauge coupling unification. The doublet-triplet splitting is achieved 
through a generalization of the missing vev mechanism. The proton decay rate is suppressed by 
arranging, without fine-tuning, a pair of Higgs doublets to lie at an intermediate scale. Gauge 
coupling unification is then restored by an appropriate splitting of the SU(5) components of 
the 54 containing the type II seesaw triplet, which in turn brings a$(mz) within la from 
its experimental value, thus improving on the MSSM prediction. In passing, we provided an 
updated, comprehensive analysis of proton decay from D=5 and D=6 operators. 
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A SO (10) gauge symmetry breaking 



The purpose of this appendix is to provide an explicit sector breaking the SO(10) gauge 
symmetry down to the SM gauge group. It is well known that this breaking can be realized by 
the most general renormalizable superpotential involving the following Higgs representations: 
one 45, one 54 and one (16, 16) pair. In this case, both SM singlets in the 45 acquire a 
nonzero vev. However, we need a 45 multiplet with a vev aligned in the B — L direction to 
implement the missing vev mechanism for doublet-triplet splitting (see Appendix [Bj. In order 
to obtain such a vacuum alignment, we must introduce additional fields and consider a non- 
generic superpotential. A simple realization of this is provided by the following superpotential: 

W S b = ^/iS^^Si + (A 12 S + /i 2 54 , )45i452 + ^A54 , 54 / 54 / +T6(M 16 + (745016 , (A.l) 

where S is an SO(10) singlet and the contractions of SO(10) indices are understood. As 
shown below, the 54' acquires a GUT-scale vev and therefore cannot be identified with the 54 
multiplet involved in the type II seesaw mechanism. 

Altogether, eight SM singlets can acquire a nonzero vev at the GUT scale. We normalize 
their vevs as follows: (Tr(54't54')) = |Vp 5 | 2 , which breaks SO(10) down to its Pati-Salam 
subgroup, (Tr(45j45i)) = \V®_ L \ 2 + \V®\ 2 (i = 1,2), (16tl6) = |Ui| 2 , (16*16) = |Fi| 2 and 
(S) = S. In order to preserve supersymmetry, these vevs must satisfy F-flatness and D-flatness 
conditions. The solution to these constraints readi§ 



(A.2) 



together with |Vi| = |Vi|. If the dimensionless couplings are of order one, then all nonzero vevs 
are of order M\q, the unique mass parameter in Wsb- Therefore SO(10) is broken down to the 
SM gauge group in one step. The two 45 vevs are aligned along the T 3 jj and B — L directions, 
respectively; in the rest of the paper we will rename them Vr = vjp and Vb-l = Vb~l- 

The SM vacuum defined by Eq. (|A.2p provides all uneaten chiral superfields with a GUT- 
scale mass, except for a pair of charged states with SU(3)c x SU (2)lxU (l)y quantum numbers 
(1, 1, ±1), which remains masslesqj. This problem can be cured by adding the following terms 
to W S b: 

- A 3 S 45 3 45 3 + Ai 23 45i45 2 45 3 . (A.3) 

This does not modify the vacuum alignment discussed above and does not introduce extra 
mass parameters either. The F-flatness conditions imply that the 45 3 does not acquires a vev. 
The coupling Ai 23 is sufficient to make the unwanted massless states heavy, while the coupling 
A 3 guarantees that all 45 3 components are also massive. 



s Let us mention for completeness that the D-flatness and F-flatness conditions admit another solution char- 
acterized by Vg}_ L = Vjp = and all other vevs nonzero. This solution would also be satisfactory for our 
purposes, but we stick to the first one for definiteness. 

9 In the alternative vacuum characterized by Vg_ L = = 0, the massless states are (3, 1, — 2/3)©(3, 1,2/3). 
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For completeness, we give the masses of the X and Y gauge bosons, which can mediate 
proton decay through D = 6 operators (see Appendix IB. 3|) : 

M\ = Ml = g 2 GVT (| \V PS \ 2 + I \V B ^ L \ 2 + \ |1/ R | 2 ) , (A.4) 

where <7gut is the SO(10) gauge coupling. 

We remark that the desired vev alignment is obtained because several couplings allowed by 
the SO(10) gauge symmetry are absent in Eq. (|A.1|) . In fact, one can add mass terms for 45i 
and 54' without upsetting the vev alignment. On the contrary, it is crucial to forbid all 452 
couplings except A12 and f\%. We do not try to justify the non-generic form of VFsb by global 
symmetries, which would require a more complicated set of fields and interactions. 

Note finally that it is possible to align the vev of one of the two 45 's without introducing 
the SO(10) singlet S. In fact, if one replaces A12 S in Eq. f| A. 1|) by a bare mass M12, the 
F-flatness equations still have a solution with = and a second one with v£l L = 0, all 
other vevs being nonzero. This is sufficient to realize the doublet-triplet splitting, but this 
requires the presence of two mass parameters (Mi 6 and M12) in Wsb- 

B Doublet-triplet splitting and proton decay 

In this appendix, we provide an explicit doublet-triplet splitting mechanism for the SO (10) 
scenario studied in this paper, and derive the conditions imposed by the non-observation of 
proton decay on the heavy spectrum. 

B.l Doublet-triplet splitting 

In the SO(10) implementation of the type II seesaw mechanism considered in this paper, up 
quark masses arise from the superpotential couplings 16^16^10, while down quark and charged 
lepton masses arise from the couplings 16jl0jl6. The Y = +1/2 MSSM Higgs doublet h u 
should therefore reside mainly in the 10 in order to account for the large top quark mass, while 
the Y = —1/2 Higgs doublet should contain a sig nificant Hf component, and we require 
that it also has a non-vanishing Hi component!^. The doublet-triplet splitting mechanism, 
besides rendering all colour triplets coupling to light matter fields heavy, should therefore leave 
and an admixture of Hh and massless. 
This can be achieved by a generalization of the missing vev mechanism [29J involving an 
additional 10' Higgs multiplet and an adjoint Higgs field with a non-vanishing vev in the B — L 
direction (which motivates the choice made for 452 in Appendix |A|) . The superpotential that 
accomplishes the desired doublet-triplet splitting reads: 

W DT = -Mio'10 / 10 / + /ilO / 45 2 10 + T6(Mi 6 +545i)16 + -T6T6(fylO + plO / ) , (B.l) 

where the third term also plays a role in SO(10) symmetry breaking (and therefore appears in 
Eq. (lA.ljl ). while the coupling p 16 16 10' is optional. As in the case of Eq. (|A7l1) . this is not 
the most general superpotential for the fields involved. After SO(10) symmetry breaking, one 

10 Indeed, the leptogenesis scenario of Ref. [7] assumes that the dominant decay modes of the heavy matter 
fields preserve B — L, which is guaranteed if hd contains a non-negligible H^ component. 
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obtains the following doublet and triplet mass matrices, written in the (10, 10', 16) (left) and 
(10, 10', T6) (right) bases: 

/ o o -Wi \ ( o ^ v b-l -Wi \ 

M D = M w -pVx , M T = ~^Vb-l M w -pV x . (B.2) 

V Mi 6 / \ Mi 6 + 2gV R ) 

From Eq. (|B.2|) we can see that all colour triplets and two pairs of Higgs doublets acquire 
GUT-scale masses, while and a combination of and remain massless: 

K = H™ , h d = cos6 H Hf + sm6 H Hf , (B.3) 

where the Higgs mixing angle Oh is given by (with no loss of generality, we assume fjV\ and 
Mi6 to be real): 

tanflff = ^ . (B.4) 

Note that the alignment of the 452 vev along the B — L direction is crucial for the splitting of 
the doublet and triplet components. 

A few comments are in order about the above doublet-triplet mechanism. First, the SO(10)- 
breaking sector contains a 54' representation with a vev in the Pati-Salam direction, whose 
couplings to 10 10 and 10 10' must be forbidden as they would make all Higgs doublets heavy. 
On the contrary, the coupling 10 10 54 involved in the seesaw mechanism is harmless since 54 
does not acquire a vev. Second, the mass term Migl616 in Eqs. (jA.ip and (jB.ip is crucial both 
for GUT symmetry breaking and for the doublet-triplet splitting: if it were absent, F-fiatness 
would imply either V\ = 0, leaving the SO (10) rank unbroken, or Vr = 0, leaving a pair of 
colour triplets massless. 

B.2 Higgs-mediated proton decay 

In general, giving GUT-scale masses to the colour Higgs triplets is not enough to suppress 
the proton decay rate below its experimental upper limit. To achieve this, one must impose 
additional constraints on the triplet mass matrix. We show below that this can be done in a 
simple way in the doublet-triplet splitting scenario presented above. 

Let us first adapt the standard computation of the Higgs-mediated proton decay rate [111 
1431 HH1 H7] to our case. Upon integrating out the heavy Higgs triplet superfields, one obtains 
the following D=5 operators: 



(QiQj)iQkLi) 



+ o K 'ijkl u i u jdk e l 
d 2 ^ 



+ h.c, (B.5) 

e 2 



where the contraction of gauge indices is understood, and the dimensionful coefficients Kjiki 
and K'j ikl generated by the superpotential (P) are given bJ^l: 

Hijki = ~ Vijhki (Mf 1 ) T10? i6 , n' ijkl = - (yuhjk ~ VjiKk) (Mf 1 )^^ , (B.6) 

where (My 1 ) Tl0 — ie is the (T 10 ,T 16 ) entry of the inverse triplet mass matrix. Since colour 
invariance implies Kmi = and K' iik , = 0, the dominant proton decay modes arising from 



11 Note that the flavour structure of Kijki and Kyjy is the same as in the minimal SU(5) model. This is an 
obvious consequence of the way the SM matter fields are embedded into SO(10) representations. 
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the above operators involve a kaon, and in practice p — > K + u dominates. The corresponding 
amplitude is obtained by "dressing" the D=5 operators of Eq. (|B.5j) with gaugino/higgsino 
loops. If the first two generation squarks are close in mass, as the observed weakness of 
hadronic flavour- violating processes may suggest, the dominant contributions come from the 
wino dressing of the QQQL operator and from the charged higgsino dressing of the u c u c d°e c 
operator. Here we consider only the former contribution; the latter (which is significant for 
the K + v T channel only) would give a stronger constraint on Mt only for large values of tan/3. 
The proton partial decay width then reads: 



r(p -» k a 
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32vr/2 m 3 
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1 + 



m r , 
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(D - F) + 



6mA 



(D + 3F) 



C\2\l 



(B.7) 



where f n = 130 MeV is the pion decay constant; (3h is the hadronic parameter defined by 
(3hPlu p = {0\(ulcLl)ul\p) , where u p is the proton spinor and the parenthesis indicates the con- 
traction of Lorentz indices; D and F are chiral Lagrangian parameters; and Cijki is the Wilson 
coefficient of the four-fermion operator (uudL^idLk^i)- The most recent lattice determination 
of Ph is \Ph\ = 0.0120(26) GeV 3 [48], while an analysis of hyperon decay measurements gives 
F + D = 1.2670 ± 0.0030 and F — D = -0.341 ± 0.016 @9]. The Wilson coefficients C ijkt read: 



Cijkl — ^ ^ VmjVnk(Kimnl 



K"nmil j 



f{d' Li ,u Lm ) + f(u Ln ,e L i) 



(B.8) 



where An is a renormalization factor, V is the CKM matrix, the k coefficients are expressed 
in the basis q = (u, d') = (u, Vd), and the loop function / is given by (M2 is the wino mass): 
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For degenerate sfermion masses (m a = m^ = in) and a hierarchy Mf <C m 2 , f(a,b) reduces 
to iV/2/m 2 to a good approximation. In Eq. (IB.8j) . the quantity (Ki mn i — K nrn u) is evalu- 
ated at the GUT scale, where the operator QQQL is generated. The renormalization factor 
Ar = AsdAld contains a short-distance piece A$d which accounts for the renormalization 
of the superpotential operator QQQL from Mqut to the supersymmetry breaking scale (here 
identified with mz), and a long-distance piece Ald which encodes the renormalization of the 
four-fermion operator (uLid£,j)(dife^) from mz to Hhad = 1 GeV. The latter is given by: 
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and the former by (neglecting the Yukawa contributions): 
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1.4, 



(B.10) 



(B.ll) 



where / runs over all mass thresholds M/ between m^ and MqtjTj and 6{ 2 3 are the beta 
function coefficients between Mj and Mj + i. Given the hierarchy among the quark Yukawa 
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couplings at the GUT scale, one has Civu — Cuii, with: 

C xm - f^i? (M r 1 ) T10Tl6 dlul l ' J2 e2inX ^VndV ns \f(d' L ,u Ln ) + f(u Ln ,e Ll )] , (B.12) 



4tt " V 1 'T^T- siQ H 

where <E>^ and <3?^ (n = 1,2,3, 3>" + ^2 + ^3 = 0) are high-energy phases, and the sum is 
dominated by the n = 2 term (i.e. by the &l loops). 

We are now in a position to derive a lower bound on (MjT ) Tl0 — 16 from the experimental 
constraint r(p — > K + v) > 2.3 x 10 33 yrs (90% C.L.) [50]. For definiteness, we consider the 
model (ii) of Section 14.31 with Ma = 10 12 GeV, tan /3 = 10 and a spectrum close to the one of 
Fig. |6(b)[ with squark masses in the TeV range (except for the lightest stop), slepton masses 
around 750 GeV and M 2 « 280 GeV. This spectrum gives f(d' L , c L ) + f{c L , e u ) « 0.45 TeV -1 . 
Evaluating the Yukawa couplings at the GUT scale, we find X u V u dV us = 5.2 x 10~ 7 , X c V c dV cs = 
-2.6 x 10~ 4 , Xt\V td V t s\ = 7.5 x 10" 5 , and {X d V ud , X S V US , X b \V ub \) = (0.43,1.9,2.0) x 10~ 4 . 
Depending on the values of the high-energy phases, a partial cancellation between the ti and 
the cl contributions is possible. Here we do not consider this possibility and keep only the 
dominant cl contribution, which gives: 

\C im \ 2 - 9.9 x lO-^GeV' 1 ( ^) ( ( ^— j) (M^Y 10 - 6 , (B.13) 

1 V 8.1 J \sme H J V0.45TfeV _ V Jt 10 t ' v J 

where we have replaced f(d' L ,c£) + /(cL,ey) by an average value 2/, and A$d — 8.1 (to be 
compared with 7.8 in the MSSM with tan/3 = 10) takes into account the various thresholds 




shown in Fig. 1(a) From Eqs. (|B.7|) and (|B.13p . we finally obtain: 



2.3xl0 33 yrs V Iftrl /V 1 J\ V J\ [M^^ie 



t( P ^K+u) / 0.012 GeV 3 V /sin ^V/0.45 TeV" 1 \ 2 / (4.4 x 10 18 GeV) 



(B.14) 



which yields the upper bound (Mf 1 )^^ < [4.4 x 10 18 GeV(l/sin0#) (tan/3/10) ] , where 
we have restored the approximate tan (5 dependence of the proton decay amplitude (one can 
indeed see from Eq. (|B.12j) that C\\n and Cvm scale as 1/ sin 2/3). This bound is a conserva- 
tive one, and several effects could relax it: there are still large uncertainties on the hadronic 
parameter /3#; A$d overestimates the running of the k coefficients, since it does not include 
the Yukawa contribution; a superpartner spectrum characterized by heavier sfermions and/or 
a stronger sfermion/gaugino mass hierarchy would reduce the proton decay amplitude; correc- 
tions to the mass relation M d = Mj could affect the couplings of the Higgs triplets to quarks 
and leptons; and finally, some sfermion mixing patterns compatible with the observed level of 
flavour violation can significantly reduce the proton decay amplitude [51]. 

Let us now translate this bound into constraints on the doublet-triplet splitting parameters. 
From Eq. (|B.2j) we have: 

, n Vi (VEfjM 10 , - phV B _ L ) 3 tan fl H M 1(y y^pFi 

lTW6 ' V2M m (hV B „ L )i {hV B - L ? V2M m hV B - L ' [ Ab) 

where we made use of Vr = M\§j2g and tan Oh = f\V\jM\§. Since the SO(10) gauge symmetry 
is broken in one step (see Appendix|A|), V\ ~ Vb-l ~ -^16 ~ Mqut and the desired suppression 
of (My 1 ) Tl0 ^i6 with respect to its natural value Mq^j T can be achieved by a partial cancellation 
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in the combination (V677.M10' — pIiVb-l), or by taking p <C 1 and Mic in the (10 13 — 10 14 ) GeV 
range, where the lower value corresponds to the conservative bound on (M^ 1 ) Tl0 ^i6. Note 
that tan#jj <C 1 would not help, since the upper bound on (My 1 )^^ becomes stricter for 
smaller values of Oh- In this paper we takd 1 " 2 ! M\$> = 10 14 GeV in order to avoid a strong 
cancellation among unrelated superpotential parameters. One can check that all colour Higgs 
triplets acquire GUT-scale masses in this case, while one pair of Higgs doublet sits at the 
intermediate scale M 10 '. The consequences for gauge coupling unification are discussed in 
Appendix [Cj and summarized in Section [4.31 



B.3 Gauge-mediated proton decay 

Let us now discuss the contribution of SO(10) gauge interactions to proton decay [17]. Due 
to the way the SM matter fields are embedded into SO(10) representations, there are no new 
gauge contributions with respect to the SU(5) case, i.e. only the X and Y gauge bosons 
mediate proton decay. The dominant decay mode from (X,Y) exchange is p — ► 7r°e + , and its 
rate is given by (neglecting the non-gauge contributions, which are subdominant): 

F(p -» A+) = {m } ~ f*f \a H \ 2 (1 + D + Ff {\C RL \ 2 + \C LR \ 2 ) , (B.16) 

where an is the hadronic parameter defined by ohPlu p = {0\(urcIr)ul\p) , D and F are the 
same chiral Lagrangian parameters as above, and Crl and Clr are the Wilson coefficients 
of the four-fermion operators (urcIr) {ul &l) and (uLd,L)(uReR), respectively. The most recent 
lattice determination of an is \otji 

I = 0.0112(25) GeV 3 [48J. The Wilson coefficients take the 

same form as in minimal SU(5): 

Crl = A R ^ , C LR = A R ^ (l + |K,| 2 ) , (B.17) 

where My = Mx = My is the mass of the heavy (X, Y) gauge bosons. The long-distance piece 
of the renormalization factor Ar = A$dAld is given by Eq. (IB.lOjl . while its short-distance 
piece reads: 

23 3 4 

i-r ( ai(Mj) \ ( q 2 (Mj) \ W 2 ( a^Mi) \ iSJ 

AsD ~ 1/ \MM^)) V^(Mm)J \MM^T)) ' ( } 

where, as in Eq. (jB.llh . the Yukawa contributions have been neglected, and I runs over all 
mass thresholds Mi between mz and Mqtjt- Together with Eq. (|B.17p . Eq. ()B.16P gives: 

o , . , S4 N / 0.0112 GeV 3 \ 2 / 2.4 \ 2 / 1/24 \ 2 / My \ 4 . 



to be compared with the experimental upper bound t( P -» vr°e + ) > 8.2 xlO 33 yrs (90% C.L.) [52 
In Eq. (|B.19p . Agp = 2.4 and ojgut = 1/24 are the reference MSSM values. Consider now 
the model (i) of Section with M A = M E = 10 12 GeV, M T = 10 13 GeV, M H = 10 14 GeV 



12 Note that the mass scale M 10 i ~ 10 14 GeV can be generated by an interaction term A SlO'lO' with A ~ 10 
where S is the SO(10) singlet introduced in Appendix lAl 
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and the same heavy fermion spectrum as in Fig. 1(a) At the one- loop level and omitting 
both low-energy and GUT thresholds, one obtains ogut = 1/14, Mqut = 3.5 x 10 15 GeV 
and A$d = 2.4. These values are at odds with the experimental limit on the proton lifetime 
(assuming My = Mqut leads to r(p — > 7r°e + ) = 4.2 x 10 32 yrs). On the contrary, for model (ii) 
with the spectrum of Fig. [1(a)] (with M A = 10 12 GeV, M s = M s = M T = M = 10 13 GeV 
and M H = 10 14 GeV), one finds ogut = 1/12, M GU T = 1.2 x 10 16 GeV and A SD = 2.5, 
leading to r(p — * 7r°e + ) = 3.9 x 10 34 yrs for My = Mqut- We shall therefore adopt model (ii) 
as our reference model, although model (i) may be viable if 2-loop running and threshold 
effects conspire to increase the GUT scale. Also, the corrections needed to depart from the 
minimal SU(5) relation M^ = will in general modify Eq. (|B.17[) by introducing non-trivial 
fermion mixing angles in the heavy gauge boson couplings, and this could significantly reduce 
the proton decay rate 



B.4 Proton decay from non-renormalizable operators 

For completeness, we mention that proton decay could also be induced by non-renormalizable 
superpotential operators of the form: 

16,16^10*16, (B.20) 

where A is the cutoff, and the SU(5)-singlet component of the 16 acquires a GUT-scale vev. 
These operators, if present, will generate the D = 5 operators of Eq. (|B.5|) after SO(10) 
symmetry breaking. In order to avoid a confict with the experimental bound on proton lifetime, 
they must either be forbidden by some symmetry, or the ones involving light generation fields 
must be suppressed by small coefficients dijki- This is actually what one would expect in a 
theory of flavour capable of explaining the hierarchy of fermion masses. Note that the required 
suppression of the coefficients is less severe than in conventional SO(10) models, where the 
dangerous operators, of the form 16,16^16^16;, have dimension 5. 



C Intermediate scales and gauge coupling unification 

In this appendix, we study the constraints imposed by the requirement of successful gauge 
coupling unification on the extra heavy states present below the GUT scale, using 1-loop 
renormalization group equations. 

As shown in Appendix lB.21 the experimental constraint on the p — > K + u rate can easily be 
satisfied by allowing a pair of Higgs doublets to lie at an intermediate scale. This tends to spoil 
unification, but we will show that the problem can be cured by splitting the components of 
the 54. In fact, it is also desirable to give GUT-scale masses to some components of the SU(5) 
multiplets (15, 15) and 24 contained in the 54 in order to preserve perturbative unification. 
Indeed, the presence of additional chiral superfields at intermediate scales increases the value 
of the unified gauge coupling ogut with respect to the MSSM. One has to check that a Landau 
pole is not reached soon above (or even below) the GUT scale, because this would make the 
predictions of the model very sensitive to the effects of higher-dimensional operators. In the 
SO(10) scenario studied in this paper, keeping the full 54 close to the seesaw scale would 
give c^uT — 1/4, and the Landau pole would be reached at a scale A 2Mgtjt- Thus, 
motivated both by perturbativity and by proton decay, we are led to split the (15, 15) and 
24 multiplets, keeping below the GUT scale only the components that are necessary for the 
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seesaw mechanism and for leptogenesis, together with possible additional components needed 
to achieve unification. 

As is customary, we define Mqtjt as the scale where a\ and «2 meet. Denoting by M n 
the masses of the intermediate-scale states and by b\ n ^ (i = 1,2,3) their contributions to the 
beta-function coefficients, the 1-loop predictions for a^{mz), Mqut and «gut read: 
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where the superscript "0" refers to MSSM quantities, with (&?, b% b^) = (^,1,-3), Mg UT 
2 x 10 16 GeV, a° GVT ~ 1/24, and 



bi - bj , 



(C.4) 



Given the fact that the MSSM prediction for a^mz) significantly deviates from the mea- 
sured value (using 2-loop RGEs and including low-energy supersymmetric thresholds, one 
finds [l/ag Xp (mz) — l/a^(mz)] — +4.3/(2-7r), which corresponds to a 5<r deviation), a positive 
contribution of the extra states to Eq. (|C.1|) would be welcome. 
Replacing Eq. (lG2l) into Eq. (ftX3l) one finds: 



«gut 



"GUT 



&2 , Mqut 
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M„ 



(C.5) 



The third term always increases ckgut> while the second one can decrease or increase it, de- 
pending on whether Mqut is smaller or larger than Mq UT . However, the contribution of the 
second term is bounded by the requirement Mqut ^5x 10 15 GeV coming from proton decay 
(see Apprendix IB.3|) . Therefore, in order to avoid a too large value of «gut> the additional 
intermediate-scale fields besides the ones needed to realize leptogenesis and to suppress pro- 
ton decay [namely (A, A), S and/or T, (5j,5j) and (H,H)\ should better be SU(2)l singlets. 
There are two such 54 components: (X, X) and O. Since O is an electroweak singlet, it does 
not affect Mqut nor «gut (but it corrects the prediction for a^). Adding only O would push 
Mqut above the Planck scale. We are thus left with two possibilities: (i) intermediate (X, X); 
(ii) intermediate (X, X) and O. The right-hand side of Eq. (|C.1|) reads, for each of the two 
cases: 
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while Mgut and «gut are given by the same expression in both cases: 
M GUT 1 



In 



M° UT 15 
1 1 



CKQUT a GUT ^'-* 7r 



_ 8 i n Mhx + ln Mhx + 5 ln Mhx + ln ^gut 
M E M A M T M H 

88 l n MUX _ n ]n ^GUT _ 85 ln ^GUT _ 2 g ln Mut 

M s M A M t Mff 

-i5m. (m - t)3 



(C.7) 



M 5l M 52 M 53 . ■ (( ' ,N) 

where we assumed M^ = Md? = Ms 4 . From the above equations, we can see that lowering 
Ms decreases qgut as desired (and improves the prediction for 03), but decreases Mgut- 

Let us first consider case (i). Assuming M E = M A = M 15 = 10 12 GeV, M T = M 24 = 
10 13 GeV and M H = 10 14 GeV (as well as \ H = 0.045, Vi = M GU T, tan j3 = 10 and tan 0# = 1 
to fix the masses of the heavy (5j, 5j) pairs), we obtain, at the 1-loop level: 

1 1 U J!) M GUT = 3.5 x 10 15 GeV , — !— = 14.1 . (C.9) 



a 3 (m z ) al(m z ) 2tt ' ' o G ut 

In this case, the prediction for a 3 (m z ) is only slightly larger than in the MSSM, and the value 
of the unified coupling remains reasonable. Unfortunately, the unification scale lies almost 
one order of magnitude below the MSSM prediction, which leads to a too fast proton decay 
rate (see Appendix IB.3j) . Gauge coupling unification is approximately preserved if M15 and 
M24 are varied while keeping the ratio M24/M15 fixed. In particular, increasing M15 while 
keeping M24/M15 = 10 slightly increases Mgut (and decreases «3 and oqut)- For instance, 
for M15 = 10 13 GeV, one obtains l/a 3 (m z )-l/a%(rn z ) = +0.27/(2vr), M GU T = 4.8xl0 15 GeV 
and ckgut = 1/15.9. Still the unification scale is dangerously close to the proton decay bound. 

Let us now consider case (ii). As anticipated, only a 3 {m z ) is affected by the presence of 
O below Mgut- This gives the possibility of increasing Mgut with respect to case (i) by 
increasing Ms, while correcting the prediction for a 3 {m z ) by adjusting Mo- For instance, the 
choice M A = 10 12 GeV, M T = M = M E = 10 13 GeV and M H = 10 14 GeV gives, at the 
1-loop level: 

T \ ~ 0/ ^ = ' M «ut = 1.2 x 10 16 GeV , = 12.5 . (CIO) 

a 3 {m z ) a^{m z ) 2tt «gut 

In this case, unification works better than in the MSSM. Indeed, we checked that the pre- 
diction for a 3 {m z ), including low-energy thresholds and the 2-loop MSSM running, matches 
the measured value within la. Moreover, there is no conflict between the value of Mgut and 
proton decay, and the Landau pole lies one order of magnitude above Mgut- From Eq. (|G6b ). 
we can see that the contributions of T and O cancel if My = Mq. Therefore, unification is 
still preserved if the masses of the various states are varied while keeping M7 1 = Mo and 
M^/ Ml » 10 15 GeV. 

The splitting of the masses of the 54 components needed to realize case (ii) can be achieved 
with the set of operators shown in TableEl where Vps is the Pati-Salam invariant vev of the 54', 
Vb-l is the vev of the 452, which is aligned along the B — L direction, and 453 has no vev (see 
Appendix|A]). At the renormalizable level, only (Z, Z) and (V, V) acquire a (GUT-scale) mass. 
The dimension-5 operators provide Ms = My = Mo ~ M GUT /A as well as Ms ~ M GUT /A 
(more precisely, when the GUT-scale masses of the 453 and 54' components are taken into 
account, Ms,t,o and Ms are further suppressed by a mild seesaw-like mechanism). Finally, 
the dimension-6 operator generates M A ~ M GUT /A 2 . 
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Operator 


Massive 54 components 


Mass 


54.453.54' 
—(54.453)45(16.16)45 

j (54.45 2 ) 5 4(54 , .54 , )54 
^2 (54.54)i(l6.45 2 .16)i 


(Z,Z),(V,V) 
0, 1 , u 

(£,£) 

(A, A) 


V PS 
-V lVl 

-^Vb-lVps 



Table 3: The operators listed in the first column generate masses for the 54 components given 
in the second column. The order of magnitude of these masses is reported in the third column. 



D Renormalization group equations 

Below MquT) the superpotential terms ([2]) and ([3]) read, in terms of the SM components 

W = (X u )ij u\qjh u + (X d )ij d c i q j h d + (A e )y efljh d 

+ (X^ij D? qj h d + (A e )y e\Ljh d + (M^ijLiLj + (M D c)^i)?L»| 



+ -L (/o)^ d?05| + (/ T ) 4J LjTZj + A {fsi)ij liSlj - -L (/sjo- %SDt 

1 /"IF" 1 1 - 

+ -j= a T h u Th d + y ^ °"5 ^S/id + -(7a ft-dA/id + — o" A /i u A/t u 

+ M A A A + M Z ZZ + M 2 SS + M V VV + -(M 5 S 2 + M T T 2 + M O 2 ) , (D.l) 

where Clebsch-Gordan coefficients have been factorized out, and contractions of SU(3)c and 
SU(2)l indices are understood. The interactions of fields with GUT-scale masses, namely 
the right-handed neutrinos and the components of 16 and 10 other than the light Higgs dou- 
blets, have been omitted^!. The boundary conditions at the GUT scale for the superpotential 
couplings and mass parameters are: 

X u = y , X e = X d = sin 9 H h , X e = X d = cos H y , 

f x = f for X = A,A,Z,Z,i:,fl,V,V,0,T,Si,S d , 

cta = cos 2 6 ho , 0"A = c , ar = &s = cosOhp , (D-2) 

M L = M D c = hVi , 

where the Higgs mixing angle 6h is defined by Eq. (|B.3|) . and V\ is the vev of the SU(5) 
singlet component of the 16. We did not write boundary conditions for the masses of the 54 



Similarly, the interactions of the 54 components with GUT-scale masses should not appear in Eq. (|D.1|I 
nor in the RGEs. 
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components, since, as discussed in Appendix O they are assumed to be split by operators not 
included in the superpotential ([3j). 

The soft supersymmetry breaking terms for the SO(10) fields are defined by: 

-L soft = (m? 6 )yl6|16j + 04)^-10*10; + m? 6H 16*16 + m 2 W[I 10*10 + m| 4 54*54 

+ Q(A 2/ )yl6 i 16 i 10 + {A^ijlSilOjie + ^(A / ) ii 10 i 10 i 54 + ^10 10 54 + h.c.^ 

+ ^M 1/2 A a A a + h.c.j , (D.3) 

where we used the same notation for the chiral multiplets and for their scalar components, A a 
are the SO(10) gauginos, and we omitted the S-terms. The soft terms of the SM components 
are given, at the GUT scale, by the following boundary conditions: 

2 2T 2T 2 2T 2 

m q = m u c = m e c = m L = m D c = m 16 , 

2 2T 2 2T 2 

m l = m d c = = m^, c = m 10 , 

m h u = m io H > m l d = cos 2 H m 2 10H + sm 2 6 H m 2 l6[i , (D.4) 
m\ = m 2 bA for X = A, A, Z, Z, S, S, S, T, O, V, V , 
Mi = M 2 = M 3 = M 1/2 . 

The boundary conditions for the ^4-terms, not included in the above list, are analogous to the 
ones for the corresponding Yukawa couplings. One has for instance: 

A u = A y , A e = A% = sm6 H A h , A e = A d = cos 6 H A y . (D.5) 

In the numerical study of Section [5l we solved the 1-loop RGEs for all superpotential 
couplings and soft terms below Mqut- For brevity, we only list below the RGEs for the 
MSSM parameters, which are sufficient for a leading-log analysis of flavour-violating effects. 
Let us first recall the 1-loop RGEs for gauge couplings and gaugino masses: 



1 - l -b a , b a = -3C 2 (G a ) + Y,T a (R) , (D.6) 



d 

dt a 2ir 

R 



^-M a = ^-b a a a M a , (D.7) 
at 2tt 

where t = log(/// / uo)> being the renormalization scale and a reference scale, C2(G a ) is the 
second Casimir invariant of the group G a , T a (R) is the Dynkin index of the representation R, 
and the sum in b a runs over all chiral superfields with mass smaller than \x. 
The 1-loop RGEs for the MSSM Yukawa couplings are: 

(4tt) 2 ^A u = X u (3\i\ u + X j d \ d + X d Xd) + Tr ^3A U A^ X u 

+ (j Wt\ 2 + ^ \as\ 2 + \ k A | 2 ) A u - C^gl + 3g% + ^g 2 ) X u , (D.8) 
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(47r) 2 -A d = X d (3X% + X[X U + 3^) + TV (3\ d \l + A e A^ + X e Xl + 3A d Aj) X d 
+ [fzf\ + fvfl + 2/e/J + \fof ] o + ^/5 d /J d ) A, 

+ (| kr| 2 + ^ ks| 2 + | kA| 2 ) A, - + 3^ + A d , (D.9) 

(4vr) 2 ^A e = (3A e At + 3A e At) A e + Tr (3A^ + A e At + A e At + 3A^) A e 
+ A e Q/I/a + |/J/z + + + 

+ (| kT| 2 + ^ks| 2 + | |a A | 2 ) A e - Q 5 2 + 3 5 2 ) A e . (D.10) 
The 1-loop RGEs for the MSSM soft sfermion masses are: 

(4vr) 2 ^m 2 = (m 2 q XiX u + AjjA M m 2 + 2X\ 1 X u m 2 Lu + 2Xim 2 uC X u + 2A{A^ 
+ (m 2 q X\X d + X\X d m 2 q + 2X\X d m 2 hd + 2X\m 2 dc X d + 2A\A d ^j 
+ (m 2 A+A rf + A* A d m 2 + 2X\X d m 2 hi + + 2i+i d ) 

- (^il V + 6|M 2 | 2 5 2 + f |M 3 | 2 5 f - \ 9 IS) , (D.ll) 

(4vr) 2 ^m 2 c = 2(m 2 c A^ + A (i AX c + 2A ( iAX d + 2A (i m 2 A^ + 2A (i 4) 
+ (m% fzfl + /z/lm 2 c + 2/ z 4m| + 2f z m 2 f\ + 2^^ z ) 
+ 2 (m 2 c /s/| + /s/|m 2 c + 2/ s 4m 2 + 2/ s m 2 ^ + 2A^A^) 
+ ^ (n&fo/J, + /5 d 4m 2 c + 2/ Sd 4m 2 + 2f Sd m%fl + 2A /Srf 4j 
+ £ + foflml + 2/o/Jm 2 , + 2f m%& + 2A fo A\ o ) 

+ ("4/v/ y + fvflml + 2/y/Jm 2 , + 2/ y m|/ y + 2A }v A) v ) 

/ \ / 8 32 2 \ 

+ 2 (m 2 cDc A,A^ + X d X d m 2 d l Dc ) - ^-|Mx| 2 5l 2 + - |M 3 | 2 ff 2 - -<7?5j , (D.12) 

(4vr) 2 ^m 2 c = 2(m 2 c A u At +A u Atm 2 c + 2A u Atm 2 u + 2A u m 2 At +2A u At) 

-g|M 1 | 2 5 2 + ||M 3 | 2 5 2 + ^) , (D.13) 
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(4vr) 2 -mf = (mfxlX e + XtX e mf + 2XlX e ml d + 2Xlmlc\ e + 2AlA e j 

+ \ (mfflU + fif A mf + 2f{f A m 2 A + 2f{m 2T f A + 2^ A/a ) 

+ I (mifzfz + + 2flf z m 2 z + 2/Jm3c/z + 2A f /z A /z ) 

+ \ (mtfvfv + 4 W + VUWv + 24m^ c /v + 2A)_ A f . ) 

+ | (W/t/t + 4W + 24/rr4 + 2/^mf / r + 2^^) 

+ ^ + + 2/4/ Si m| + 2/J jTO f f SL + 2A} s 

+ (At X e mf L + m 2 L X\ A e ) - (j\MA 2 g 2 + 6|M 2 | 2 ff 2 + jj 5 2 ^ , (D.14) 



(4vr) 2 -m 2 c = 2(m 2 c A e At + A e Atm 2 c + 2A e Atm 2 d + 2A e m 2 At + 2^t) 
+ 2 (m 2 c A e At + A e A];m 2 c + 2X e X\m\ d + 2X e m 2 L Xl + 2A e A^ 

-(^\Mi?!ii-\<ils) , (D.15) 

where the hypercharge D-term contribution S is given by: 

S = m\ u - m\ d + Tr(-m 2 + m 2 c - 2m 2 c + m% + m 2 ) + Tr(-m| + m 2 ^ + m\ - m\ c ) 
+ 3(t7& - m\) + (m| - m|) + 4(m| - m|) + 5(m 2 / - m^) . (D.16) 

The mixing terms m 2 d c Dc d c * D c +m 2 L l* L+\\.c. appearing in the RGEs for m 2 dc and mf, although 
absent at Mgut, are generated by the RGE evolution and must be included for consistency. 

Finally, the 1-loop RGEs for the MSSM A-terms are: 

(4tt) 2 ^U u = A u (bXlXu + X\X d + X%) + 4A u At A u + 2A U (x\A d + \\A d ) 
+ Tr (3A u At ) A u + \a T \ 2 + ^ \a s \ 2 + | |<r A | 2 ) A u 

+ Tr (QA U X\^ X u + {^A aT a* T + j^^s^S + 3At a ^a) A « 

+ — g\ (2M 1 X U - A u ) + 3 52 2 (2M 2 A M - A u ) + — g 2 (2M 3 X U - A u ) , (D.17) 
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(47r) 2 -^ = A d (5\l\ d + At \ u + 5Al A d ) + A\ d \ ] d A d + 2\ d (x[A u + 2X ] d A d ) 
+ [fzfl + fvfl + 2/e/J + pofl + ^fs d fl d ) A d 
+ + 2A fv fl + 4A fs fl + S -A f J Q + ^A fs / S )j \ d 

+ Tr (3A d A d + A e At + A e A| + 3A d A d ) A d + Q |<r T | 2 + ^ |a 5 | 2 + | |<t a | 2 ) 4* 
+ 2Tr (3A*Aj + A e At + i e At + 3i d A rf ) A rf + QAtt^t + ^ A °s°*s + ^A aA a^ X d 
+ ^g\ (2M 1 X d - A d ) + 3g 2 2 (2M 2 X d - A d ) + ™$ (2M 3 X d - A d ) , (D.18) 

{^?j t Ae = 4 (A e Al + A e At) A e + 5A e X\ X e + 5A e X\ X e 

+ A e Q/t / A + |/t /z + |/t /y + 3/t fT + A/t ; / 5 ^ 

+ Ae (sflA fA + 3/| A /z + 3/t A fv + \&A fT + ^/|// S( ) 

+ Tr (3A d A d + A e At + A e At + 3A d A d ) A e + Q |a T | 2 + A | CT5 | 2 + 1 j^ 2 ^ ,4 e 

+ 2Tr (3A d A d + A e A* + i e A* + 3i d A d ) A e + {^-A aT a* T + ^A as a* s + 3A aA <j* A ^ X e 

+ \g\{2M 1 X e - A e )+3g 2 2 {2M 2 X e - A e ) . (D.19) 
5 
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